Method for producing complex real three-dimensional images, and system for same

ABSTRACT

The present invention relates to a method for producing complex reality three-dimensional images and a system for same, the method comprising: (a) a step for determining first reality three-dimensional spatial coordinates for a three-dimensional image of a human body; (b) a step for determining second reality three-dimensional spatial coordinates for an image of an item of medical equipment; (c) a step for obtaining a three-dimensional image of the area surrounding the medical equipment, from an imaging means in the medical equipment, and determining third reality three-dimensional spatial coordinates for said image; (d) a step for examining an image that is at the same coordinates in the three kinds of three-dimensional spatial coordinates; and (e) a step for producing a complex reality three-dimensional image by selecting the one image that is at the same coordinates, if there is one image at the same coordinates, or selecting the necessary image or images from among a plurality of images, if there are multiple images at the same coordinates.

TECHNICAL FIELD

The present invention relates to a method for producing a complex reality three-dimensional image and a system for the same.

BACKGROUND ART

During operations in which a medical device is invasively or non-invasively inserted into the human body, a precise spatial location of the medical device inserted into the human body as well as precise spatial locations of human tissues is required in order to increase the success rate of surgery and prevent unexpected medical emergencies.

In the conventional medical technology, for the diagnosis and treatment of a patient, a medical device that enters the body of the patient is inserted without a precise location thereof after three dimensional images of the body of the patient are captured in advance using magnetic resonance imaging (MRI), computed tomography (CT), PET CT, gamma camera, or the like.

This medical technology has problems in that tissues of the patient body may be unnecessarily dissected and the undesired tissue dissection may cause fetal complications since three-dimensional spatial coordinates of the patient body are known through MRI or CT but three-dimensional coordinates of the medical device inserted into the patient body cannot be grasped in real time during operations.

In order to solve these problems, there is provided “a relative coordinate extracting device and a medical imaging device using the same (Korean Patent Publication No. 10-2012-0096729)” as a technology of obtaining spatial coordinates of the patient body during operations. This technology discloses that an image capturing device is installed in the patient body and relative coordinates of the image capturing device with respect to a target marker are calculated, but has problems in that accurate spatial coordinates of the patient cannot be recognized in a region of the human body in which the image capturing device is installed, and in cases where a prompt surgical operation needs to be done due to the critical condition of the patient, the installation of the image capturing device and the coordinate analysis may cause waste of so called golden time.

In addition, an example of a technology for analyzing a three-dimensional image of a particular subject is “a three-dimensional shape measurement method and a three dimensional measurement system (Korean Patent Publication No: 10-2012-0090964”. This technology discloses that coordinate values representing a three-dimensional shape of a first part of a subject to be measured and coordinate values representing a three-dimensional shape of a second part of the subject are converted into coordinate values based upon a reference coordinate system identified by marks, and the coordinate values after the conversion are synthesized, so that entire three-dimensional shape data of the subject are obtained. This technology is somewhat similar to a proposal of the present invention with respect to three-dimensional spatial coordination of a medical device. However, this technology requires a separate light source and a reference coordination system, which cannot be used at the time of separate operations, and thus has a decisive problem in that this equipment cannot be placed in an operation room. Moreover, this technology merely recognizes the three-dimensional shape of a particular subject, but cannot describe the coordinates of the particular subject in a three-dimensional space. Therefore, this technology cannot be applied to the three-dimensional space coordinates of the medical device. Recently, it is a trend to perform operations and procedures using various images, and here, the excellence of performance cannot be provided since the operations or procedures are performed while different images are determined through different imaging devices.

Throughout the entire specification, many cited documents and patent documents are referenced and their citations are represented. The disclosures of cited documents and patent documents are entirely incorporated by reference into the present specification, and the level of the technical field within which the present invention falls, and details of the present invention are explained more clearly.

DETAILED DESCRIPTION OF THE INVENTION Technical Problem

The present inventors have endeavored to solve the problems of the conventional art. As a result, the present inventors have developed a novel method capable of accurately finding a three-dimensional position of a medical device used in a procedure, by allowing three-dimensional spatial coordinates of the human body and three-dimensional spatial coordinates of the medical device to coincide with each other to integrate different kinds of images of the human body. According to the present invention, a reality three-dimensional image capable of grasping an actual procedure in real time in a procedure place can be produced.

Therefore, the present invention has been made in view of the above-mentioned problems, and an aspect of the present invention is to provide a method for producing a complex reality three-dimensional image.

Another aspect of the present invention is to provide a system for producing a complex reality three-dimensional image.

Other purposes and advantages of the present invention will become clarified by the following examples, claims, and drawings.

Technical Solution

In accordance with an aspect of the present invention, there is provided a method for producing a complex reality three-dimensional image, the method including:

(a) determining first reality three-dimensional spatial coordinates of a three-dimensional image of the human body obtained before a medical device is applied to a subject;

(b) determining second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject, the medical device including an imaging unit;

(c) obtaining a three-dimensional image of the periphery of the medical device from the imaging unit of the medical device, and determining third reality three-dimensional spatial coordinates of the obtained image, the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates being determined by using the same coordinate system;

(d) examining images on the same coordinates in the three three-dimensional spatial coordinates by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates; and

(e) producing a complex reality three-dimensional image by selecting, if one image is on the same coordinates, the one image, and selecting, if multiple images are on the same coordinates, necessary image(s) from the multiple images.

The present inventors have endeavored to solve the problems of the conventional art. As a result, the present inventors have developed a novel method capable of accurately grasping a three-dimensional position of a medical device used in a procedure, by allowing three-dimensional spatial coordinates of the human body and three-dimensional spatial coordinates of the medical device to coincide with each other and integrating different kinds of images of the human body. According to the present invention, a reality three-dimensional image capable of grasping an actual procedure in real time in a procedure place can be produced.

The present invention is directed to a complex reality three-dimensional image, and provides a feature of producing a single three-dimensional image by incorporating different images.

The images incorporated herein are a three-dimensional image of the human body, a three-dimensional image of the medical device, and a three-dimensional image of the periphery of the medical device (e.g., an image obtained by the image medical device). Preferably, the images incorporated herein are a three-dimensional image of the human body, a three-dimensional image of the medical device, and a three-dimensional image of the periphery of the medical device, and an image linked with signals of the human body or an image obtained by an image medical device. As used herein to refer to the three-dimensional image, the term “complex” refers to the incorporation of the different images into one image. As used herein to refer to the three-dimensional image, the term “reality” refers to showing an actual situation, in which a procedure is performed, in real time through the produced image. As used herein to refer to the three-dimensional spatial coordinates, the term “reality” refers to three-dimensional coordinates in a space or environment in which a procedure is performed.

The present invention has been described by the method for producing a complex reality three-dimensional image, but the present invention may be described by a method for unifying spatial coordinates of a three-dimensional image of the human body and spatial coordinates of a three-dimensional image of the medical device.

According to the present invention, in a space other than a space in which a procedure (e.g., an operation) is performed, spatial coordinates of a three-dimensional image obtained by scanning the human body through an imaging equipment are secured before a medical device is applied to a subject, and in the space in which the procedure is performed, three-dimensional spatial coordinates of a particular part of the human body are calculated by using various imaging methods, for example, X-ray fluoroscopy. The calculated three-dimensional spatial coordinates of the particular part are allowed to coincide with the coordinates of the particular part on the previously secured three-dimensional spatial coordinates of the human body, thereby obtaining first reality three-dimensional spatial coordinates of the three-dimensional image of the human body, which is a subject of the procedure. As for the medical device, a two-dimensional image of the medical device is secured by various image methods, for example, X-ray fluoroscopy, and then, second reality three-dimensional spatial coordinates of the medical device (including, e.g., a medical unit for working and a medical imaging unit) based on the distance relationship among the medical device, an X-ray light source, and an X-ray detector. Third three-dimensional spatial coordinates of a three-dimensional image obtained by the medical imaging unit are determined using the three-dimensional spatial coordinates of the medical device. Then, the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body, the second reality three-dimensional spatial coordinates of the image of the medical device, and the third reality three-dimensional spatial coordinates of the three-dimensional image obtained by the medical imaging unit are incorporated to reconstitute a necessary image. Thus, the present invention can provide accuracy and convenience of the procedure.

According to an embodiment of the present invention, the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates may be on an orthogonal coordinate system.

According to an embodiment of the present invention, the same coordinate system used to determine the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates may be set based on an operation table on which the medical device is applied to the subject.

According to an embodiment of the present invention, the three-dimensional image of the human body used in step (a) may be an X-ray fluoroscopy image, a computed tomography (CT) image, a magnetic resonance imaging (MRI) image, a positron emission tomography (PET) image, a PET/CT image, a PET/MRI image, a radioisotope imaging (RI) image, an ultrasonography image, or a complex image thereof.

According to an embodiment of the present invention, the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates may be determined using electromagnetic radiation. The electromagnetic radiation is X-ray radiation or radiofrequency electromagnetic radiation, and particularly X-ray radiation.

According to an embodiment of the present invention, the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates may be determined using X-ray fluoroscopy.

According to an embodiment of the present invention, the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates may be determined using mono-plan fluoroscopy, bi-plan fluoroscopy, or multi-plan fluoroscopy.

According to an embodiment of the present invention, the first reality three-dimensional spatial coordinates of step (a) may be determined using three-dimensional spatial coordinates of the three-dimensional image of the human body obtained before the medical device is applied to the subject and reality three-dimensional spatial coordinates of a particular part of the subject human body.

According to an embodiment of the present invention, step (a) may include the following sub-steps of: (a-1) irradiating X-rays toward the particular part of the human body of the subject using a first X-ray light source; (a-2) detecting an X-ray fluoroscopy image of the human body by the first X-ray light source, using a first X-ray detector; (a-3) irradiating X-rays toward the particular part of the human body to cross the X-rays by the first X-ray light source, using a second X-ray light source; (a-4) detecting an X-ray fluoroscopy image of the human body by the second X-ray light source, using a second X-ray detector; and (a-5) determining the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body by using the human body images detected by the first X-ray detector and the second X-ray detector.

According to an embodiment of the present invention, sub-step (a-5) may include: (a-5-1) determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector, by using physical three-dimensional coordinate information of the first X-ray light source and the first X-ray detector, distance information therebetween, and size information of the particular part of the human body detected by the first X-ray detector, and determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector, by using physical three-dimensional coordinate information of the second X-ray light source and the second X-ray detector, distance information therebetween, and size information of the particular portion of the human body detected by the second X-ray detector; (a-5-2) determining three-dimensional spatial coordinates of the particular part of the human body, by using the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector and the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector; and (a-5-3) allowing the three-dimensional spatial coordinates of the particular part of the human body to coincide with the three-dimensional spatial coordinates of the three-dimensional image of the human body of step (a) to determine the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body.

According to an embodiment of the present invention, the medical device may include an identification part for recognizing the position of the medical device. Since the entire three-dimensional shape of the medical device shows a relative position based on the identification part, a three-dimensional image of the medical device is implemented by using the identification part in a step of implementing a complex reality image.

According to an embodiment of the present invention, the identification part may include at least one of: (i) at least one sphere-shaped identification part, (ii) at least one ring-shaped identification part, (iii) an identification part including a ring and at least one marker indicating a specific site of the ring, (iv) an identification part including at least one pair of rings, (v) an identification part including at least one pair of rings and at least one marker indicating a specific site of the rings, or (vi) an identification part including markers respectively indicating three differentiated points.

Particularly, the three points of (vi) the identification part are differentiated by using markers with different images, indicating the respective points, or indicating the points using ends of linear markers having differentiated images or intersection points of the markers.

According to an embodiment of the present invention, the identification part may include at least one pair of rings and at least one marker. Here, the pair of rings are mounted on the medical device at different positions, and the at least one marker is attached to a position (one side of the ring) at which the position of the ring is capable of being indicated.

According to an embodiment of the present invention, the at least one pair of rings may have different imaging characteristics. The term “different imaging characteristics” refers to some characteristics (e.g., position, thickness, shape, and radiation transmittance) that make images be differentiately shown. For example, the rings with different thicknesses are differentiately shown through X-ray fluoroscopy.

According to an embodiment of the present invention, in cases where the number of markers is two or more, the markers may have different imaging characteristics (e.g., position, thickness, shape, length, and radiation transmittance). According to a more specific embodiment, the two or more makers have different lengths, and thus different images are calculated to give different imaging characteristics.

According to an embodiment of the present invention, the second reality three-dimensional spatial coordinates may be determined using one or two or more X-ray fluoroscopy images, and step (b) may include the following sub-steps of: (b-1) irradiating X-rays toward the medical device using an X-ray light source; (b-2) detecting an X-ray fluoroscopy image of the medical device by the X-ray light source, by using an X-ray detector; and (b-3) determining the second reality three-dimensional spatial coordinates of the image of the medical device, by using the through-position of the identification part, which is seen through on the X-ray fluoroscopy image.

According to an embodiment of the present invention, in sub-step (b-3), the second reality three-dimensional spatial coordinates of the image of the medical device and three-dimensional spatial coordinates of the marker may be determined by using the through-position of the identification part, which is seen through on the X-ray fluoroscopy image, thereby grasping the relative position of an internal portion of the medical device from the three-dimensional spatial coordinates of the marker.

According to an embodiment of the present invention, sub-step (b-3) may be performed using three-dimensional position information of the identification part, or three-dimensional position information and dimension information (e.g., diameters of rings, distance between rings, and a shape of the ring attached on the medical device). More specifically, in sub-step (b-3), the second reality three-dimensional spatial coordinates of the image of the medical device are determined by using physical three-dimensional coordinate information of an X-ray light source and an X-ray detector, distance information therebetween, size information of the medical device detected by the X-ray detector, a diameter of the ring, and the distance between the rings.

According to an embodiment of the present invention, the medical device may include a radio wave reflector. According to an embodiment of the present invention, step (b) may include the following sub-steps of: (b-1) generating a first microwave toward the radio wave reflector using a first microwave generator; (b-2) generating a second microwave toward the radio wave reflector using a second microwave generator; (b-3), after sub-steps (b-1) and (b-2), receiving the first and second microwaves reflected from the radio wave reflector; and (b-4) determining the second reality three-dimensional spatial coordinates of the image of the medical device by using information of the received first and second microwaves.

According to an embodiment of the present invention, the imaging unit included in the medical device may include a plurality of cameras, and here, the three-dimensional image of the periphery of the medical device is obtained using the plurality of cameras in step (c). According to a more specific embodiment, the camera as the imaging unit may include an identification part, and more specifically, a pair of rings and a marker attached to one side of each of the rings, respectively, as the identification part.

According to an embodiment of the present invention, the imaging unit may be an ultrasonic probe, and here, the three-dimensional image of the periphery of the medical device is obtained using the ultrasonic probe in step (c). According to a more specific embodiment, the ultrasonic probe as the imaging unit may include an identification part, and more specifically, a pair of rings, a marker attached to one side of each of the rings, and a scan line, as the identification part.

According to an embodiment of the present invention, the three-dimensional image of the human body of step (a) and the three-dimensional image of the periphery of the medical device of step (c) may be different kinds of images.

According to an embodiment of the present invention, the image of the medical device and the three-dimensional image of the periphery of the medical device, which are on the same coordinates and selected in step (e) may be images at the same time point.

According to an embodiment of the present invention, the method of the present invention may further include a step of determining a three-dimensional image of a periodic-mobility human body part (e.g., heart), which has periodic mobility and generates signals (e.g., electrocardiogram, blood pressure, organ pressure, and abdominal pressure) linked with the periodic mobility, and fourth reality spatial coordinates thereof, and here, step (d) may be performed by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, the third reality three-dimensional spatial coordinates, and the fourth reality spatial coordinates.

In accordance with another aspect of the present invention, there is provided

A system for producing a complex reality three-dimensional image, the system including:

(a) a first process for calculating first reality three-dimensional spatial coordinates of a three-dimensional image of the human body obtained before a medical device is applied to a subject; (ii) second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject; and (iii) third reality three-dimensional spatial coordinates of a three-dimensional image of the periphery of the medical device, by using the same coordinate system; and

(b) a second processor for examining images on the same coordinates in the three three-dimensional spatial coordinates, by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial, and the third reality three-dimensional spatial, and then, selecting, if one image is on the same coordinates, the one image, and selecting, if a plurality of images are on the same coordinates, necessary image(s) from the plurality of images, to produce a complex reality three-dimensional image.

The system of the present invention is for implementing the foregoing method of the present invention, and descriptions of overlapping contents therebetween are omitted to avoid excessive complication of the present specification.

According to an embodiment of the present invention, the system may further include an X-ray fluoroscopy imaging device, a computed tomography (CT) imaging device, a magnetic resonance imaging (MRI) imaging device, a positron emission tomography (PET) imaging device, a PET/CT imaging device, a PET/MRI imaging device, a radioisotope imaging (RI) imaging device, or an ultrasonography imaging device, for obtaining the three-dimensional image of the human body.

According to an embodiment of the present invention, the system further include an X-ray fluoroscopy imaging device or a radiofrequency electromagnetic radiation device for obtaining the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates.

According to an embodiment of the present invention, the system may further include a medical device, and the medical device may include an identification part for recognizing the position of the medical device.

According to an embodiment of the present invention, the medical device used in the present invention may include an imaging unit.

Advantageous Effects

Features and advantages of the present invention are summarized as follows:

(a) According to the present invention, an image capable of showing a three-dimensional movement of the medical device in real time in a three-dimensional image of the human body, which is obtained by allowing spatial coordinates of the human body in a place, where a medical practice is performed, to coincide with coordinates of a previously secured three-dimensional image of the human body and a three-dimensional image obtained from an imaging medical device on the place, can be produced.

(b) According to the present invention, the space of the medical device inserted into the human body can be recognized by grasping and coordinating the three-dimensional space of the medical device such that the previously secured three-dimensional image of the human body is allowed to coincide with a space location of the human body in the place where a medical practice is performed.

(c) According to the present invention, the three-dimensional space of the medical device can be recognized in real time through X-ray scans employed in the operation place.

(d) A precise operation can be performed through the medical device by recognizing three-dimensional coordinates of the medical device and the human body and then incorporating these coordinates.

(e) Since three-dimensional spatial coordinates of the medical device and the human body can be accurately recognized using X-ray light sources, microwaves, ultrasonography images, and 3D cameras, unnecessary dissection of internal organs of the patient and unexpected medical accidents can be prevented based on the accurate location of the human body and the medical device.

(f) Since the spatial coordinates of the medical device can be recognized through an X-ray scanner, which is necessarily used in the operation place, unnecessary external equipments need not be delivered into the operation place.

(g) Since the medical device operates inside the human body by allowing accurate three-dimensional spatial coordinates of the medical device to coincide with three-dimensional coordinates of the human body, a more precise and safer operation can be performed.

(h) Since the three-dimensional spatial coordinates of the human body, the three-dimensional spatial coordinates of the medical device, and the three-dimensional spatial coordinates of images of the human body and the medical device captured by an imaging medical device are combined to display the three-dimensional image, a more precise operation can be performed.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a view showing that X-rays are irradiated to a subject (a medical device) using a plurality of X-ray light sources.

FIG. 2( a) shows the position of point A projected on a detector and the position of the detector in a three-dimensional space where three axes of an orthogonal coordinate system are plotted, and FIG. 2( b) shows point A on the detector on a two-dimensional plane of a rectangular coordinate system.

FIG. 3 is a view showing the plurality of X-ray light sources and an image of the subject in FIG. 1, and the projection of the subject on a three-dimensional plane.

FIG. 4 is a conceptual perspective view showing a cylinder-shaped medical device attached with an identification part (rings and a marker).

FIG. 5 is a view showing an image of a medical device on an X-ray detector and the projected medical device when X-rays are irradiated to the medical device from one X-ray light source P in a three-dimensional space (here, an identification part is omitted).

FIG. 6 is a three-dimensional plane view showing a portion (represented by hatching) of the medical device, through which X-rays pass from the X-ray light source P, wherein the portion forms an image of one ring and an image of straight line AB connecting edges of the ring image on the detector, in FIG. 5.

FIG. 7 is a view showing the X-ray light source P, the portion (represented by hatching) of the medical device through which X-rays pass, and straight line AB (virtual line) of the image on the detector, on a two-dimensional plane. The X-ray light source P, the portion of the medical device, and straight line AB are present on one plane in FIG. 6.

FIG. 8 is a view showing a medical device and an image of the medical device, when it is assumed that edge points of images of two rings are respectively A, B, C, and D, an intersection point of the straight line bisecting ∠APB and straight line AB is E, and an intersection point of the straight line bisecting ∠CPD and straight line CD is F.

FIG. 9 is a view showing points G′, H′, I′, and J′ on a three-dimensional plane. In FIG. 8, when a straight light connecting points E and F is extended and the intersection points of the straight line and the images of the respective rings are G, H, I, and J, theses points are formed by passing the X-rays projected from the X-ray light source P through the center axis of the medical device including the rings. Here, points on the rings through which the X-ray pass are G′, H′, I′, and J′.

FIG. 10 shows a relationship among the points G, H, I, and J, which are formed by passing the X-rays through the center axis of the medical device, points G′, H′, I′, and J′ corresponding to these, and point P on a two-dimensional plane. Points G, H, I, and J, points G′, H′, I′, and J′, and point P are present on the same plane in FIG. 9. The medium point connecting points G′ and H′ is represented by O_(a), and the medium point connecting points I′ and J′ is represented by O_(b).

FIG. 11 shows a medical device including a marker indicating one point of a ring, an X-ray light source P, and an image of the medical device on a detector. The origin point of the ring is represented by O_(a); the intersection point of the ring and the marker is represented by M′; and the image thereof is represented by M. Points G, H, I, and J of FIG. 10 are shown.

FIG. 12 is a perspective view of an external shape of a camera attached with a pair of rings and bar-shaped markers indicating a particular site of the ring.

FIG. 13 is a view showing images of a subject focused on 3D cameras.

FIG. 14 is a view showing monitor images of a subject and images of the subject focused on image sensors.

FIG. 15 is a view showing schematic straight lines of the image in FIG. 12 using three-dimensional coordinates.

FIG. 16 is a schematic view showing an ultrasonography imaging medical device with an identification including rings, markers, and a scan line.

MODE FOR CARRYING OUT THE INVENTION

The method and system for producing a complex reality three-dimensional image by incorporating spatial coordinates of the human body and a medical device, according to the present invention, can be variously modified, and may have several embodiments, and thus particular examples will be illustrated on the drawings, and detailed descriptions will be set forth. However, it shall be noted that it is not intended to limit the present invention to specific exemplary embodiments but intended to cover all the modifications, equivalents, or substitutions belonging to the technical idea and technical scope of the present invention.

Hereinafter, the method for producing a complex reality three-dimensional image of the present invention will be described with reference to drawings.

Step (a): Determining First Reality Three-Dimensional Spatial Coordinates of Three-Dimensional Image of Human Body

According to the present invention, first reality three-dimensional spatial coordinates of a three-dimensional image of the human body, which is obtained before a medical device is applied to a subject, are determined.

The method for producing a complex reality three-dimensional image according to the present invention is carried out by obtaining first reality three-dimensional spatial coordinates of a three-dimensional image of the body of a subject (specifically, a human being, more specifically, a human patient) and second reality three-dimensional coordinates of an image of a medical device, respectively, and then allowing these three-dimensional spatial coordinates to coincide with each other.

The three-dimensional image of the human body obtained before the medical device is applied to the subject, which is used in the present invention, is obtained through various imaging equipments. The human body is scanned by an appropriate imaging equipment, thereby producing a three-dimensional image. The three-dimensional image of the human body, certainly, includes images of internal organs and bones as well as the skin.

According to an embodiment of the present invention, the three-dimensional image of the human body is an X-ray fluoroscopy image, a computed tomography (CT) image, a magnetic resonance imaging (MRI) image, a positron emission tomography (PET) image, a PET/CT image, a PET/MRI image, a radioisotope imaging (RI) image, an ultrasonography image, or a complex image thereof. According to a specific embodiment, a CT image or an MRI image is preferable.

The three-dimensional image of the human body may be one image obtained from the imaging equipment. Selectively, the three-dimensional image of the human body may be a complex image made by combining multiple images.

The first reality three-dimensional spatial coordinates may be determined using various electromagnetic rays, specifically, X-rays or high-frequency electromagnetic rays, and more specifically, X-rays.

In cases where the first reality three-dimensional spatial coordinates are determined using X-rays, the first reality three-dimensional spatial coordinates may be determined using X-ray fluoroscopy (mono-plan fluoroscopy, bi-plan fluoroscopy, or multi-plan fluoroscopy).

The first reality three-dimensional spatial coordinates of the three-dimensional image of the human body in step(a) may be determined by using the three-dimensional coordinates of the three-dimensional image of the human body obtained before the medical device is applied to the subject and reality three-dimensional spatial coordinates of a particular part of the subject human body. Alternatively, the three-dimensional spatial coordinates of the three-dimensional image of the human body may be directly determined without using the reality three-dimensional spatial coordinates of the particular part of the human body. For example, in the current procedure practice, an equipment for obtaining the three-dimensional image of the human body is installed on an operation table, and the three-dimensional image of the human body is obtained directly from a patient present on the operation table, thereby determining the reality three-dimensional spatial coordinates of the three-dimensional image of the human body. Alternatively, the reality three-dimensional spatial coordinates may be determined using the three-dimensional image obtained from the operation table and a three-dimensional image obtained from outside the operation table.

The reason is that, it is very effective to calculate and use reality three-dimensional spatial coordinates of a three-dimensional image obtained from CT scans or the like using X-ray fluoroscopy since the three-dimensional image of the human body from CT or MRI scans is obtained using a contrast agent, and thus has a different image quality from an image obtained from X-ray fluoroscopy without a contrast agent, and the contrast agent cannot be continuously used during an operating procedure.

The X-ray fluoroscopy may be performed in various manners. Specifically, the determining of the first reality three-dimensional spatial coordinates includes the following steps of: (a-1) irradiating X-rays toward the particular part of the human body of the subject using a first X-ray light source; (a-2) detecting an X-ray fluoroscopy image of the human body by the first X-ray light source, using a first X-ray detector; (a-3) irradiating X-rays toward the particular part of the human body to cross the X-rays by the first X-ray light source, using a second X-ray light source; (a-4) detecting an X-ray fluoroscopy image of the human body by the second X-ray light source, using a second X-ray detector; and (a-5) determining first reality three-dimensional spatial coordinates of the three-dimensional image of the human body by using the human body images detected by the first X-ray detector and the second X-ray detector.

Here, sub-step (a-5) includes: (a-5-1) determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector, by using physical three-dimensional coordinate information of the first X-ray light source and the first X-ray detector, distance information therebetween, and size information of the particular part of the human body detected by the first X-ray detector, and determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector, by using physical three-dimensional coordinate information of the second X-ray light source and the second X-ray detector, distance information therebetween, and size information of the particular portion of the human body detected by the second X-ray detector; (a-5-2) determining three-dimensional spatial coordinates of the particular part of the human body, by using the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector and the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector; and (a-5-3) allowing the three-dimensional spatial coordinates of the particular part of the human body to coincide with the three-dimensional spatial coordinates of the three-dimensional image of the human body of step (a) to determine the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body.

The first reality three-dimensional spatial coordinates of the particular part of the human body, which can be well identified by X-ray fluoroscopy, are calculated by using a plurality of X-ray fluoroscopies, and reality three-dimensional spatial coordinates of the other parts are calculated based on the first reality three-dimensional spatial coordinates. The reality three-dimensional spatial coordinates calculated by using X-ray fluoroscopy can have a decreased error in cases where several parts of the human body are calculated. If the three-dimensional spatial coordinates of the three-dimensional image obtained by CT or MRI scans in an operating space are identical to the reality three-dimensional spatial coordinates, this process may be omitted.

The first reality three-dimensional spatial coordinates of the three-dimensional image of the human body, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates may be determined by using various coordinate systems capable of expressing three dimensions. For example, an orthogonal coordinate system may be used, but another coordinate system, for example, a spherical coordinate system or the like may be used as a reference coordinate system according to the environment in which the present invention is implemented, used equipments, and the like.

The same coordinate system is used in order to allow the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates to coincide with one another, and in this case, various points of a space in which the procedure is performed may be used as references. Among the points, it is appropriate to use the operation table as a reference in the aspect of convenience.

For example, the coordinate system using the operation table as a reference may be constructed as follows. The operation table is horizontal, and the side at which the head of a patient is placed is an upper side, and the side of the leg of the patient is a lower end. When the left and right sides are set based on the observation from the lower end to the upper end, the x-, y-, and z-axes and origin point O are set as follows: The x-axis is a line which is parallel with the short axis of the operation table and horizontal while passing through the lower end of the operation table; the y-axis is a line which is vertical and meets the x-axis while passing through a point at which the X-ray detector or the X-ray light source is positioned at the left side as much as possible; and the z-axis, which is parallel with the long axis of the operation table and horizontal while passing through an intersection point of the x-axis and the y-axis.

According to an embodiment of the present invention, in cases where the operation table is used as the reference, the gradations on the operation table, which are helpful to set three-dimensional coordinates, are favorable in setting three-dimensional coordinates.

Step (b): Determining Second Reality Three-Dimensional Spatial Coordinates of Image of Medical Device

Next, second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject are determined. The medical device includes an imaging unit.

Apart from three-dimensional imaging of the human body and spatial coordinating of the three-dimensional image, the medical device is spatially coordinated.

Here, the medical device includes all medical devices which are in contact with an external part of the patient body or inserted into the human body. The medical devices used herein may be classified into two kinds, and one is an imaging medical device and the other is a functioning medical device. The functioning medical device is a tool for performing a procedure, such as scissors, a knife, and a needle, including a stent. Specifically, the medical device used herein includes various medical devices, such as an endoscope, an ultrasonography imaging device, and a surgical robot.

The caution here is that there is no order in time between the three-dimensional spatial coordination of the human body and the three-dimensional spatial coordination of the position of the medical device. Therefore, claims indicate the spatial coordination of the human body in step (a) and the spatial coordination of the medical device in step (b), but the spatial coordination of the medical device may be first performed, and this shall be construed to be included within the scope of a right of the present invention.

The three-dimensional spatial coordination of the human body and the three-dimensional spatial coordination of the medical device are performed in parallel and independently. However, in general, the operation or the like is performed after CT or MRI scans of the patient, and thus in practical, the three-dimensional imaging and the three-dimensional spatial coordination of the human body precede the three-dimensional spatial coordination of the medical device.

The spatial coordinates used herein need to be used in the same coordinate system, and the units thereof need to be the same as each other. The three-dimensional spatial coordinates of the medical device may be used in, for example, an orthogonal coordinate system.

The reality three-dimensional spatial coordinates of the human body and the medical device are obtained by, for example, X-ray fluoroscopy, and an ultrasonography imaging device or a 3D camera is additionally used to assist the spatial coordination by X-ray fluoroscopy.

In cases where the spatial coordinates are obtained by X-ray fluoroscopy, X-rays are continuously irradiated from a light source during the surgery of the patient. Therefore, the X-ray fluoroscopy makes it unnecessary to use separate medical devices, and thus is very economical.

The CT scanned image is not an image of the human body which varies in real time, and thus in order to supplement this, necessary images are obtained by using and combining 3D camera or ultrasonography images. The reference point for combining images may be based on the reality three-dimensional spatial coordinates at which the human body of the patient is positioned. The 3D camera or ultrasonography image is a three-dimensional image configured based on the reference point of each device, and thus when the reality three-dimensional spatial coordinates of the reference point are known, the three-dimensional spatial coordinates of the 3D camera or ultrasonography image of the human body or the medical device can be replaced with the reality three-dimensional spatial coordinates. The reality three-dimensional spatial coordinates of the medical device including an imaging device can be calculated from the three-dimensional spatial coordinates of the reference point of the imaging device by using X-ray fluoroscopy.

In cases where three-dimensional spatial coordinates are incorporated or allowed to coincide with each other, a plurality of various image coordinates are present in the identical spatial coordinates. Of these, the coordinates of necessary images are selected to reconstitute a three-dimensional image.

For example, in cases where a stent is inserted for the treatment of aortic aneurysm, a three-dimensional image of the human body obtained from CT scans using a contrast agent is used for the shape of the aorta, and a three-dimensional image calculated by X-ray fluoroscopy is used for the medical device including the stent. Here, the medical device including the stent is positioned in the aorta, and thus the coordinates of the medical device is included in the coordinates of the aorta, so that the coordinates of the medical device are expressed by the medical device, and the coordinates of the other aorta is expressed by the aorta.

As another example, in cases where the lung surgery is performed using a three-dimensional endoscope, the air gets out of the lung at the time of surgery, and thus the operation is performed while the volume of the lung is reduced, but the three-dimensional image from CT scans does not reflect the reduced volume. Therefore, the image from a 3D camera is used for the image of the lung with a reduced volume, and the three-dimensional image from CT is used for the image of a part that is not shown through the 3D camera.

As still another example, in cases where the heart valve surgery is performed using a catheter, the catheter arrives at the heat through the blood vessel, and here, the three-dimensional image of the catheter using X-ray fluoroscopy is used for the shape of the blood vessel in the three-dimensional image from CT. The movement of the heart valve is difficult to observe by X-ray fluoroscopy, and thus an ultrasonography image (two-dimensional or three-dimensional image) is used therefor. The detailed image of the catheter is difficult to grasp through the ultrasonography image, and thus a three-dimensional image using X-ray fluoroscopy is used therefor.

In addition, second reality three-dimensional spatial coordinates of the image of the medical device applied to the subject are determined.

The determination of the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body and the second reality three-dimensional spatial coordinates of the image of the medical device will be described in detail by exemplifying X-ray fluoroscopy.

First, X-rays are irradiated toward a particular part of the human body (or a medical device) using a first X-ray light source 1 a. Thereafter, a first X-ray detector detects the X-rays that are irradiated from the first X-ray light source and passes through the particular part of the human body (or a medical device). In the same manner, X-rays are irradiated toward the particular part of the human body (or the medical device) using a second X-ray light source 1 b, and then a second X-ray detector detects the X-rays which are irradiated from the second X-ray light source and passes through the particular part of the human body (or the medical device). In this case, as shown in FIG. 1, the X-rays irradiated from the first X-ray light source and the second X-ray light source detect the same particular part of the human body (or the medical device), and cross each other.

The physical coordinate position of the first X-ray light source is information that is physically given in advance, and the physical coordinate position of the first X-ray detector is also information that is physically given in advance. In addition, the physical coordinate position of the second X-ray light source is information that is physically given in advance, and the physical coordinate position of the second X-ray detector is also information that is physically given in advance. If the particular part of the human body (or the medical device) is extremely small, the position thereof may be represented as a point at which lines meet each other, as shown in FIG. 3. In order to represent the detected particular part as a point, the particular part is a small region which is well detected by the X-ray, or the positional change of a region represented by the center point of the image needs to be small according to the detection angle. Here, the particular part of the human body is not limited to only a natural region of the human body, but includes an artificial material or natural material disposed inside or outside the human body.

In order to obtain the point at which lines meet each other in FIG. 3, real coordinates of the image on the detector needs to be obtained. When the X-ray detector is positioned in any space in the reality three-dimensional space, image A on the detector, which is generated by the X-rays, is shown in FIG. 2( a), and two-dimensional coordinates of the image on the detector are shown in FIG. 2( b).

The reality three-dimensional spatial coordinates of the position of the X-ray light source and the detector can be physically obtained using an X-ray fluoroscope installed on the operation table. Based on these, the reality three-dimensional spatial coordinates of the image on the detector may be obtained.

In FIG. 2, the center point D of the X-ray detector and the X_(d) axis and the Y_(d) axis constituting two-dimensional coordinates of the detector can be physically recognized. When the coordinates of point D are (d₁, d₂, d₃), the vector of X_(d) axis is u′=(d, e, f), the vector of Yd axis is v′=(g, h, i), and coordinates of point A on the image generated on the detector are (a, b), three-dimensional coordinates (a₁, a₂, a₃) of point A can be calculated, and the method therefor is as follows.

Point A is positioned on the straight line which passes through point H and has vector v, and distanced from point H by an interval of b. Here, point H is distanced from point D by an interval of a, and positioned on the X_(d) axis, and vector v is the same as the vector of Y_(d).

If point A is (h₁, h₂, h₃),

a=√{square root over ((d ₁ −h ₁)²+(d ₂ −h ₂)²+(d ₃ −h ₃)²)}{square root over ((d ₁ −h ₁)²+(d ₂ −h ₂)²+(d ₃ −h ₃)²)}{square root over ((d ₁ −h ₁)²+(d ₂ −h ₂)²+(d ₃ −h ₃)²)}  Equation (1)

Here, straight line X _(d) passes through point D and has

, and thus

the equation of X _(d) is given by

$\frac{\left( {x - d_{1}} \right)}{d} = {\frac{\left( {y - d_{2}} \right)}{e} = {\frac{\left( {z - d_{3}} \right)}{f}.}}$

Since point H is on straight line X _(d), the put of the coordinates of point H is as follows:

$\frac{\left( {h_{1} - d_{1}} \right)}{d} = {\frac{\left( {h_{2} - d_{2}} \right)}{e} = {\frac{\left( {h_{3} - d_{3}} \right)}{f}.}}$

Therefore,

$h_{2} = {{\frac{e}{d}\left( {h_{1} - d_{1}} \right)} + d_{2}}$ $h_{3} = {{\frac{f}{d}\left( {h_{1} - d_{1}} \right)} + {d_{3}.}}$

When these are put into equation (1),

$a = \sqrt{\left( {h_{1} - d_{1}} \right)^{2} + \left\{ {{\frac{e}{d}\left( {h_{1} - d_{1}} \right)} + d_{2} - d_{2}} \right\}^{2} + \left\{ {{\frac{f}{d}\left( {h_{1} - d_{1}} \right)} + d_{3} - d_{3}} \right\}^{2}}$ $a = \sqrt{\left( {h_{1} - d_{1}} \right)^{2} + \left\{ {\frac{e}{d}\left( {h_{1} - d_{1}} \right)} \right\}^{2} + \left\{ {\frac{f}{d}\left( {h_{1} - d_{1}} \right)} \right\}^{2}}$ $a = \sqrt{\left( {h_{1} - d_{1}} \right)^{2} + {\frac{e^{2}}{d^{2\;}}\left( {h_{1} - d_{1}} \right)^{2}} + {\frac{f^{2}}{d^{2\;}}\left( {h_{1} - d_{1}} \right)^{2\;}}}$ $a = \sqrt{{\frac{d^{2}}{d^{2}}\left( {h_{1} - d_{1}} \right)^{2}} + {\frac{e^{2}}{d^{2\;}}\left( {h_{1} - d_{1}} \right)^{2}} + {\frac{f^{2}}{d^{2}}\left( {h_{1} - d_{1}} \right)^{2}}}$ $a = {\left( {h_{1} - d_{1}} \right)\sqrt{\frac{\left( {d^{2} + e^{2} + f^{2}} \right)}{d^{2}}}}$ ${h_{1} - d_{1}} = \frac{ad}{\sqrt{d^{2} + e^{2} + f^{2}}}$ $h_{1} = {d_{1} + {\frac{d \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}.}}$

Solve it in the same manner, then

${h_{2} = {d_{2} + \frac{e \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}}},{h_{3} = {d_{3} + {\frac{f \times a}{\sqrt{d^{2} + e^{2} + f^{2\;}}}.}}}$

Thus, coordinates of point H (h₁, h₂, h₃) are given by:

$\left( {{d_{1} + \frac{d \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}},{d_{2} + \frac{e \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}},{d_{3} + \frac{f \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}}} \right).$

Assume that this point is H (h₁, h₂, h₃),

The straight line which passes through point H and has the vector of Y_(d), that is, v′(g, h, i), is given by:

$\frac{\left( {x - h_{1}} \right)}{g} = {\frac{\left( {y - h_{2}} \right)}{h} = {\frac{\left( {z - h_{3}} \right)}{i}.}}$

Since point A (a₁, a₂, a₃) is on this straight line, the following equation is given:

$\begin{matrix} {\frac{\left( {a_{1} - h_{1}} \right)}{g} = {\frac{\left( {a_{2} - h_{2}} \right)}{h} = {\frac{\left( {a_{3} - h_{3}} \right)}{i}->.}}} & {{Equation}\mspace{14mu} (2)} \end{matrix}$

Since the distance between point H and point A is b, the following equation is given:

b=√{square root over ((a ₁ −h ₁)²+(a ₂ −h ₂)²+(a ₃ −h ₃)²)}{square root over ((a ₁ −h ₁)²+(a ₂ −h ₂)²+(a ₃ −h ₃)²)}{square root over ((a ₁ −h ₁)²+(a ₂ −h ₂)²+(a ₃ −h ₃)²)}  Equation (3).

Equations (2) and (3) are summarized as follows:

$a_{1} = {h_{1} + \frac{g \times b}{\sqrt{g^{2} + h^{2} + i^{2}}}}$ $a_{2} = {h_{2} + \frac{h \times b}{\sqrt{g^{2} + h^{2} + i^{2}}}}$ $a_{3} = {h_{3} + {\frac{i \times b}{\sqrt{g^{2} + h^{2} + i^{2}}}.}}$

Then, put

${h_{1} = {d_{1} + \frac{d \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}}},{h_{2} = {d_{2} + \frac{e \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}}},{h_{3} = {d_{3} + \frac{f \times a}{\sqrt{d^{2} + e^{2} + f^{2}}}}}$

into the above equations:

$a_{1} = {d_{1} + \frac{d \times a}{\sqrt{d^{2} + e^{2} + f^{2}}} + \frac{g \times b}{\sqrt{g^{2} + h^{2} + i^{2}}}}$ $a_{2} = {d_{2} + \frac{e \times a}{\sqrt{d^{2} + e^{2} + f^{2}}} + \frac{h \times b}{\sqrt{g^{2} + h^{2} + i^{2\;}}}}$ $a_{3} = {d_{3} + \frac{f \times a}{\sqrt{d^{2} + e^{2} + f^{2}}} + {\frac{i \times b}{\sqrt{g^{2} + h^{2} + i^{2}}}.}}$

Therefore, when the center point D (d₁, d₂, d₃) of the radiation detector, vectors of reference lines of the detector, that is, X_(d) and Y_(d),

,

, and any point A (a, b) on the detector are known, three-dimensional coordinates of point A can be obtained.

After the reality three-dimensional spatial coordinates of any image on the detector are obtained as described above, coordinates of the intersection point in the three-dimensional space are obtained as shown in FIG. 3, and thus the coordinates will be coordinates of the particular part of the human body (or the medical device).

The coordinates are obtained as follows.

In FIG. 3, assuming that the image generated by shooting a subject M with an extremely small area from an X-ray light source A is represented by B and the image corresponding to a light source C is represented by D, for the sake of convenience, the equation of straight line AB through two points A and B is given by

$\frac{x - a_{1}}{a_{1} - b_{1}} = {\frac{y - a_{2}}{a_{2} - b_{2}} = {\frac{z - a_{3}}{a_{3} - b_{3}}.}}$

Assume that α₁=a₁−b₁, α₂=a₂−b₂, α₃=a₃−b₃, then

$\frac{x - a_{1}}{\alpha_{1}} = {\frac{y - a_{2}}{\alpha_{2\;}} = {\frac{z - a_{3}}{\alpha_{3\;}}.}}$

The equation of straight line CD through two points C and D is given by

${\frac{x - c_{1}}{\beta_{1}} = {\frac{y - c_{2}}{\beta_{2\;}} = \frac{z - c_{3}}{\beta_{3}}}},$

where β₁=c₁−d₁, β₂=c₂−d₂, and β₃=c₃−d₃.

Assume that the intersection point of straight lines AB and CD is M (m₁, m₂, m₃). Since point M is on straight lines AB and CD, the following equations are given:

$\begin{matrix} {\frac{m_{1} - a_{1}}{\alpha_{1}} = {\frac{m_{2} - a - 2}{\alpha_{2}} = {\frac{m_{3} - a_{3}}{\alpha_{3}}->}}} & {{Equation}\mspace{14mu} (4)} \\ {{{{\frac{m_{1} - c_{1}}{\beta_{1}} = {\frac{m_{2} - c_{2}}{\beta_{2}} = {\frac{m_{3} - c_{3}}{\beta_{3}}->{{From}\mspace{14mu} {equation}\mspace{14mu} (4)}}}},{m_{2} = {{\frac{\alpha_{2}}{\alpha_{1}}\left( {m_{1} - a_{1}} \right)} + a_{2}}}}{{{From}\mspace{14mu} {equation}\mspace{14mu} (5)},{m_{2} = {{\frac{\beta_{2}}{\beta_{1}}\left( {m_{1} - c_{1}} \right)} + c_{2}}}}{{\frac{\alpha_{2}}{\alpha_{1}}\left( {m_{1} - a_{1}} \right)} + a_{2}} = {{{\frac{\beta_{2}}{\beta_{1}}\left( {m_{1} - c_{1}} \right)} + {{c_{2}\left( {\frac{\alpha_{2}}{\alpha_{1}} - \frac{\beta_{2}}{\beta_{1}}} \right)}m_{1}}} = {{\frac{\alpha_{2}}{\alpha_{1}} \times a_{1}} - {\frac{\beta_{2}}{\beta_{1}}c_{1}} + c_{2} - a_{2}}}}{m_{1} = {\frac{{\frac{\alpha_{2}}{\alpha_{1}} \times a_{1}} - {\frac{\beta_{2}}{\beta_{1}}c_{1}} + c_{2} - a_{2}}{\frac{\alpha_{2}}{\alpha_{1}} - \frac{\beta_{2}}{\beta_{1}}}.}}} & {{Equation}\mspace{14mu} (5)} \end{matrix}$

In the same manner, m₂ and m₃ can be obtained.

As described above, radioscopy or a plurality of radioscopies used to obtain the reality three-dimensional coordinates of the particular part of the human body, and the obtained reality three-dimensional coordinates are allowed to coincide with three-dimensional coordinates of the three-dimensional image of CT or the like to obtain the reality three-dimensional spatial coordinates of the three-dimensional image of CT or the like.

Two or more fluoroscopies may be performed on the same particular part of the human body, and the more the fluoroscopies are performed, the smaller the error in the calculation of reality three-dimensional spatial coordinates.

In addition, as shown in FIG. 1, due to characteristics of X-rays spreading in all directions, the image of the particular part of the human body (or the medical device) on the detector will be enlarged as the particular part of the human body (or the medical device) is closer to the light source, and in contrast, the size of the image of the particular part of the human body (or the medical device) will be equal to an actual size of the particular part of the human body (or the medical device) as the particular part of the human body (or the medical device) is farther from the light source. The size of the image of the particular part of the human body (or the medical device) depending on the detected angle can be obtained from the three-dimensional image of CT or the like, and thus it can be calculated at which coordinates the particular part of the human body (or the medical device) is positioned between the first X-ray light source and the first detector and between the second X-ray and the second detector, by measuring the size of the image of the particular part of the human body (or the medical device). This can supplement the foregoing method in which the intersection point of two straight lines is calculated in the space to obtain the reality three-dimensional coordinates.

In the method of calculating the reality three-dimensional spatial coordinates of the medical device, the three-dimensional spatial coordinates of the medical device can be calculated more effectively when the mechanical shape and size of the medical device are known. In addition, the reality three-dimensional spatial coordinates of the medical device can also be calculated by merely a single time of X-ray fluoroscopy.

It is effective to mount an identification part on a predetermined part of the medical device and use the identification part rather than use the X-ray fluoroscopy while the entire shape and size of the medical device are known. The reason is that, in cases where the three-dimensional position of the identification part mounted on the medical device is found, the three-dimensional position of the entire medical device can be obtained by using the three-dimensional position of the identification part.

The identification part may have any shape, form, and structure as long as it is easily distinguishable from the human body, which is seen through by radiation fluoroscopy.

More particularly, the identification part is (i) at least one sphere-shaped identification part, (ii) at least one ring-shaped identification part, (iii) an identification part including a ring and at least one marker indicating a particular site of the ring, (iv) an identification part including at least one pair of rings, (v) an identification part including at least one pair of rings and at least one marker indicating a particular site of the rings, or (vi) an identification part including markers respectively indicating three differentiated points.

Of these identification parts, in a case of (vi) the identification part including markers respectively indicating three differentiated points, the three differentiated points are differentiated from one another by using makers with different images indicating the respective points or indicating the differentiated points using ends of linear markers with differentiated markers or the intersection points of the markers and the differentiated points.

According to an embodiment of the present invention, the identification part is composed of rings and markers, and the rings constitute at least one pair (the thicknesses of the rings are equal or different), and surround the medical device.

The two or more rings are spaced apart from each other on the medical device, and may be disposed at any position. For example, one ring may be positioned at one end of the medical device and the other ring may be positioned to be spaced apart from the end.

The identification may be formed of any material that has different radiation transmittance from the human body. For example, the identification part may be formed of a material having a high radiation transmittance.

According to an embodiment of the present invention, the rings may be positioned on the plane vertical with respect to the center axis of the medical device, or a plurality of markers may be mounted as needed. A plurality of markers may be mounted on the medical device, but ten or less markers are favorably disposed.

According to an embodiment of the present invention, the identification part includes two or more rings and one or more markers. The two or more rings are mounted on the medical device at different positions, and the one or more markers are mounted at positions at which the positions of the rings can be indicated. The marker is configured to indicate the specific site of the ring, and for example, the marker may be attached on one side of the ring, or the marker may indicate a specific portion of the ring even when the marker is disposed on the extension line of the linear marker without being attached to the ring at a distance.

Two or more rings may be disposed as needed, and preferably, the rings have different thicknesses so as to grasp the direction of the medical device. In addition, in cases where a marker capable of indicating the specific site of the ring is used, the relative position of an internal portion of the medical device can be grasped. The medical device may have a flexure portion, and here, in cases where two rings are disposed before and after the flexure portion, a three-dimensional flexible image can be obtained.

FIG. 4 is a conceptual perspective view showing a cylinder-shaped medical device on which an identification part (rings and a marker) is attached. Here, for the sake of convenience of calculation, two rings are positioned on a surface of a cylinder which has the same center axis with the two rings, and on a plane vertical to the center axis, and the diameters of the two rings are equal. The following drawings are the same as above.

FIG. 5 shows an image obtained by passing radiation through a medical device attached with an identification part, which is disposed at any position, on a detector, which is disposed at any position, and the medical device, using a three-dimensional plane. The position of the image is on the detector. When the radiation is irradiated from an X-ray light source P, the radiation passes through the medical device to form an image on a surface of the detector as shown in FIG. 5. The radiation irradiated from the light source P passes through contact points A′, B′, C′, and D′ of two rings to form images A, B, C, and D corresponding to the contact points A′, B′, C′, and D′ on the surface of the detector. Therefore, points A′, B′, C′, and D′ are positioned on straight lines, PA, PB, PC, and PD. Hence, lines PA, PB, and AB are disposed on the same plane. As shown in FIG. 6, this plane is a face of the radiation which obliquely cuts and passes through the cylinder-shaped medical device from the X-ray light source P. In addition, FIG. 7 shows the foregoing plane through a two-dimensional view. As shown in FIG. 7, the cut shape of the cylinder is an ellipse, and the X-ray light source P is positioned on one axis of the ellipse having center point O. Straight lines PA and PB are tangential lines of the ellipse, which pass through point P, and pass through contact points A′ and B′. In addition, assuming that the intersection point of straight line PO and straight line AB is point E, ∠APE=∠BPE. Here, since coordinates of points P, A, and B are known, the intersection point of the straight line, which passes through point P and bisects ∠APB, and straight line AB can be calculated.

Assume that ∠APE=∠BPE=α, then area S₁ of triangle ΔAPE is given by S₁=½ sin α PA× PE, and area S₂ of ΔBPE is given by S₂=½ sin α PB× PE. Thus, S₁:S₂=½ sin α PA× PE:½ sin α PB× PE. Therefore, S₁:S₂= PA: PB. Meanwhile, ΔAPE and ΔBPE have the same height and different bottom sides, S₁:S₂= AE: EB. Thus, AE: EB= PA: PB.

Assuming that the intersection points of line segments AB and CD and the lines which bisect ∠APB and ∠CPD are represented by points E and F, respectively, line segment EF is an image of the medical device obtained by passing the radiation irradiated from the X-ray light source P through the center axis of the cylinder-shaped medical device. The extension line of EF may also be the same as described above. As shown in FIGS. 8 to 10, assuming that the intersection points of the extension line of line segment EF and the images of one pair of rings are points G, H, I, and J, respectively, the radiation irradiated from the X-ray light source P passes through points G′, H′, I′, and J′ on the rings to form points G, H, I, and J on the image of the detector. Therefore, as shown in FIG. 10, points G′, H′, I′, and J′ are intersection points of a face including the center axis of the medical device, the X-ray light source P, and straight line GJ and one pair of rings mounted on the medical device.

Thus, points P, G, H, I, J, G′, H′, I′, and J′ are present on the same plane. FIG. 10 shows a two-dimensional view of this plane.

In FIG. 10, point P is an X-ray light source, and points G, H, I, and J are points of an image which is generated on the detector by passing the radiation irradiated from point P through points G′, H′, I′, and J′.

Square G′H′I′J′ is a longitudinal section of the cylinder, which is obtained by cutting the cylinder-shaped medical device between rings so as to include the center axis of the cylinder-shaped medical device. Therefore, square G′H′I′J′ is a square shape. In addition, points G′ and H′ are on the ring and face each other based on the center of on the ring. In addition, points I′ and J′ are on the ring and face each other based on the center of the ring.

Therefore, the lengths of line segments G′H′ and I′J′ are the diameters of the rings, and the distance between line segments G′I′ and H′J′ is the distance between the rings.

Meanwhile, in order to find three-dimensional spatial coordinates of the cylinder-shaped medical device, three-dimensional spatial coordinates of the ring need to be obtained. In addition, when the center axis of the cylinder is obtained based on the three-dimensional spatial coordinates of the cylinder-shaped medical device, the excellence in the operation of the medical device can be secured. Assuming that the points that bisect G′H′ and I′J′ are represented by O_(a) and O_(b), the straight line connecting the two points is the center of the cylinder-shaped medical device (not shown in FIG. 10).

Therefore, when three-dimensional spatial coordinates of G′, H′, I′, and J′ are known, the three-dimensional spatial coordinates of the medical device may be obtained.

Since point P is the X-ray light source, and points G, H, I, and J are points of an image generated on the detector by passing the radiation irradiated from point P through points G′, H′, I′, and J′, three-dimensional spatial coordinates of points P, G, H, I, and J may be obtained by using the foregoing method.

Points G′, H′, I′, and J′ are positioned on straight lines PG, PH, PI, and PJ, and thus, when the distances thereof from point P, that is, the lengths of PG′, PH′, PI′, and PJ′ are known, the three-dimensional spatial coordinates of G′, H′, I′, and J′ may be calculated by using the equation of a straight line and the equation between two points.

As shown in FIG. 10, assume that intersection points of extension lines of line segments G′I′ and H′J′ and the normal line perpendicular to the extension lines from point P are designed as S and T, respectively; the length of G′H′ and the length of I′J′ are designated as R, the length of G′I′ and the length of H′J′ are designated as L, the length of PS is m, and the length of I′S and the length of J′T are n. Here, R is the diameter of the rings, and L is the distance between the rings.

These are expressed by:

PS=m, I′S= J′T=n, ∠PSG=∠PTH′=π/2

G′H′= I′J′=R, G′I′= H′J′=L

Therefore,

PG′=√{square root over ((L+n)² +m ²)}

PH′=√{square root over ((L+n)²+(R+m)²)}{square root over ((L+n)²+(R+m)²)}

PI′ =√{square root over (n ² +m ²)}

PJ′ =√{square root over (n ²+(R+m)²)}.

Meanwhile, L and R are the distance between the rings and the diameter of the ring, and thus are previously known constants. Therefore, when m and n are found, the lengths of PG′, PH′, PI′, and PJ′ may be calculated.

Here, m and n are obtained as follows:

Assume that ∠GPH=h, ∠HPI=i, ∠IPJ=j, ∠JPS=s.

Then,

${\tan \; s} = \frac{n}{\left( {m + R} \right)}$ ${\tan \left( {s + j} \right)} = \frac{n}{m}$ ${\tan \left( {s + j + i} \right)} = \frac{\left( {L + n} \right)}{\left( {R + m} \right)}$ ${\tan \left( {s + i + j + h} \right)} = {\frac{\left( {L + n} \right)}{m}.}$

Summarize it using the formula of

${\tan \left( {\alpha + \beta} \right)} = {\frac{\left( {{\tan \; \alpha} - {\tan \; \beta}} \right)}{\left( {1 + {\tan \; \alpha \times \tan \; \beta}} \right)}.}$

Then,

$\begin{matrix} \begin{matrix} {{\tan \; j} = {\tan \left\{ {\left( {s + j} \right) - s} \right\}}} \\ {= \frac{{\tan \left( {s + j} \right)} - {\tan \; s}}{1 + {\tan \; s \times {\tan \left( {s + j} \right)}}}} \\ {= \frac{\frac{n}{m} - \frac{n}{m + R}}{1 + {\frac{n}{m} \times \frac{n}{m + R}}}} \\ {= \frac{\frac{{n\left( {m + R} \right)} - {m \times n}}{m\left( {m + R} \right)}}{\frac{{m\left( {m + R} \right)} + n^{2}}{m\left( {m + R} \right)}}} \\ {= {\frac{{mn} + {nR} - {mn}}{{m\left( {m + R} \right)} + n^{2}}.}} \end{matrix} & {{Equation}\mspace{14mu} (6)} \\ {{{\tan \; j} = \frac{nR}{m^{2} + {mR} + n^{2}}}{{m^{2} + {mR} + n^{2}} = \frac{nR}{\tan \; j}}{{m^{2} + {mR} + n^{2} - \frac{nR}{\tan \; j}} = {0->}}} & \; \\ {{{\tan \; i} = {\tan \left( {\left( {s + j + i} \right) - \left( {s + j} \right)} \right)}}\begin{matrix} {{\tan \; i} = \frac{\frac{L + n}{R + m} - \frac{n}{m}}{1 + {\frac{L + m}{R + m} \times \frac{n}{m}}}} \\ {= \frac{\frac{{m\left( {L + n} \right)} - {n\left( {R + m} \right)}}{\left( {R + m} \right)m}}{\frac{{\left( {R + m} \right)m} + {\left( {L + n} \right)n}}{\left( {R + m} \right)m}}} \\ {= \frac{{m\; L} + {mn} - {nR} - {mn}}{{\left( {R + m} \right)m} + {\left( {L + n} \right)n}}} \\ {= \frac{{m\; L} - {nR}}{m^{2} + {mR} + {nL} + n^{2}}} \end{matrix}} & {{Equation}\mspace{14mu} (7)} \\ {{{m^{2} + {mR} + {nL} + n^{2}} = \frac{{m\; L} - {nR}}{\tan \; i}}{{m^{2} + {mR} + {n\; L} + n^{2} - \frac{{m\; L} - {nR}}{\tan \; i}} = {\theta->}}} & \; \\ \begin{matrix} {{\tan \; h}\; = {\tan \left\{ {\left( {s + j + i + h} \right) - \left( {s + j + i} \right)} \right\}}} \\ {= \frac{\frac{L + n}{m} - \frac{L + n}{R + m}}{1 + {\frac{L + n}{m} \times \frac{L + n}{R + m}}}} \\ {= \frac{\frac{{\left( {L + n} \right)\left( {R + m} \right)} - {m\left( {L + n} \right)}}{m\left( {R + m} \right)}}{\frac{{m\left( {R + m} \right)} + \left( {L + n} \right)^{2}}{m\left( {R + m} \right)}}} \\ {= \frac{{\left( {L + n} \right)\left( {R + m} \right)} - {m\left( {L + n} \right)}}{{m\left( {R + m} \right)} + \left( {L + n} \right)^{2}}} \\ {= \frac{{LR} + {m\; L} + {nR} + {mn} - L - {mn}}{m^{2} + {mR} + L^{2} + {2{nL}} + n^{2}}} \end{matrix} & {{Equation}\mspace{14mu} (8)} \\ {{{\tan \; h} = \frac{{LR} + {nR}}{m^{2} + {mR} + L^{2} + {2{nL}} + n^{2}}}{{m^{2} + {mR} + L^{2} + {2{nL}} + n^{2}} = \frac{{LR} + {nR}}{\tan \; h}}{{m^{2} + {mR} + L^{2} + {2{nL}} + n^{2} - \frac{{LR} + {nR}}{\tan \; h}} = {0->}}} & \; \\ {{\tan \left( {j + i} \right)} = {{\tan \left\{ {\left( {s + j + i} \right) - s} \right\}} = {\frac{\frac{L + n}{R + m} - \frac{n}{R + m}}{1 + {\frac{L + n}{R + m} \times \frac{n}{R + m}}} = \frac{\frac{\left( {R + m} \right)L}{\left( {R + m} \right)^{2}}}{\frac{\left( {R + m} \right)^{2} + {nL} + n^{2}}{\left( {R + m} \right)^{2}}}}}} & {{Equation}\mspace{14mu} (9)} \\ {{{\tan \left( {j + i} \right)} = \frac{{RL} + {m\; L}}{R^{2} + {2{mR}} + m^{2} + {nL} + n^{2}}}{{R^{2} + 2 + m^{2} + {nL} + n^{2}} = \frac{{RL} + {Lm}}{\tan \left( {j + i} \right)}}{{R^{2} + 2 + m^{2} + {n\; L} + n^{2} - \frac{{RL} + {Lm}}{\tan \left( {j + i} \right)}} = {0->}}} & \; \\ \begin{matrix} {{\tan \left( {j + i + h} \right)} = {\tan \left\{ {\left( {s + j + i + h} \right) - s} \right\}}} \\ {= \frac{\frac{L + n}{m} - \frac{n}{R + m}}{1 + {\frac{L + n}{m} \times \frac{n}{R + m}}}} \\ {= \frac{\frac{{\left( {R + m} \right)\left( {L + n} \right)} - {mn}}{m\left( {R + m} \right)}}{\frac{{m\left( {R + m} \right)} + {nL} + n^{2}}{m\left( {R + m} \right)}}} \\ {= \frac{{RL} + {nR} + {m\; L} + {mn} - {mn}}{{m\left( {R + m} \right)} + {n\; L} + n^{2}}} \end{matrix} & {{Equation}\mspace{14mu} (10)} \\ {{{\tan \left( {j + i + h} \right)} = \frac{{RL} + {nR} + {m\; L}}{m^{2} + {mR} + {n\; L} + n^{2}}}{{m^{2} + {mR} + {nL} + n^{2}} = \frac{{RL} + {nR} + {m\; L}}{\tan \left( {j + i + h} \right)}}{{m^{2} + {mR} + {nL} + n^{2} - \frac{{RL} + {nR} + {m\; L}}{\tan \left( {j + i + h} \right)}} = {0->}}} & \; \\ \begin{matrix} {{\tan \left( {i + h} \right)} = {\tan \left\{ {\left( {s + j + i + h} \right) - \left( {s + j} \right)} \right\}}} \\ {= \frac{\frac{L + n}{m} - \frac{n}{m}}{1 + {\frac{L + n}{m} \times \frac{n}{m}}}} \\ {= \frac{\frac{{m\; L} + {mn} - {mn}}{m^{2}}}{\frac{m^{2} + {nL} + n^{2}}{m^{2}}}} \end{matrix} & {{Equation}\mspace{14mu} (11)} \\ {{{\tan \left( {i + h} \right)} = \frac{m\; L}{m^{2} + {n\; L} + n^{2}}}{{m^{2} + {n\; L} + n^{2}} = \frac{m\; L}{\tan \left( {i + h} \right)}}{{m^{2} + {n\; L} + n^{2} - \frac{m\; L}{\tan \left( {i + h} \right)}} = {0->}}} & \; \end{matrix}$

Here, equation (6)-equation (7) is given by:

$\begin{matrix} {{{{m^{2} + {mR} + n^{2} - \frac{nR}{\tan \; j} - \left( {m^{2} + {mR} + {n\; L} + n^{2} - \frac{{m\; L} - {nR}}{\tan \; i}} \right)} = {{0 - \left( {{nL} + \frac{nR}{\tan \; j}} \right) + \frac{{m\; L} - {nR}}{\tan \; i}} = 0}}{{\frac{m\; L}{\tan \; i} - {\left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right)n}} = 0}{\frac{L}{\tan \; i}m} = {\left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right)n}}{m = {{\left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right)n \times \frac{\tan \; i}{L}n}->}}} & {{Equation}\mspace{14mu} (12)} \end{matrix}$

Here, equation (8)-equation (10) is given by:

${m^{2} + {mR} + L^{2} + {2n\; L} + n^{2} - \frac{R\left( {L + n} \right)}{\tan \; h} - \left\{ {m^{2} + {mR} + {nL} + n^{2} - \frac{{RL} + {nR} + {m\; L}}{\tan \left( {j + i + h} \right)}} \right\}} = 0$ ${{nL} + L^{2} - \frac{R\left( {L + n} \right)}{\tan \; h} + \frac{{RL} + {nR} + {m\; L}}{\tan \left( {j + i + h} \right)}} = 0$ ${{nL} + L^{2} - \frac{RL}{\tan \; h} - \frac{nR}{\tan \; h} + \frac{RL}{\tan \left( {j + i + j} \right)} + \frac{nR}{\tan \left( {j + i + h} \right)} + \frac{m\; L}{\tan \left( {j + i + h} \right)}} = 0$ ${{\left( {L - \frac{R}{\tan \; h} + \frac{R}{\tan \left( {j + i + h} \right)}} \right)n} + \frac{m\; L}{\tan \left( {j + i + h} \right)} + L^{2} + {{RL}\left( {\frac{1}{\tan \left( {j + i + h} \right)} - \frac{1}{\tan \; h}} \right)}} = 0$ ${{Put}\mspace{14mu} {equation}\mspace{14mu} (12)\mspace{14mu} {in}\mspace{14mu} {{this}.\mspace{14mu} {Then}}},{{{\left( {L - \frac{R}{\tan \; h} + \frac{R}{\tan \left( {j + i + h} \right)}} \right)n} + \frac{L \times \left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right) \times \frac{\tan \; i}{L} \times n}{\tan \left( {j + i + h} \right)} + L^{2} + {{RL}\left( {\frac{1}{\tan \left( {j + i + h} \right)} - \frac{1}{\tan \; h}} \right)}} = 0}$ ${{\left( {L - \frac{R}{\tan \; h} + \frac{R}{\tan \left( {j + i + h} \right)}} \right)n} + \frac{\left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right) \times \tan \; i \times n}{\tan \left( {j + i + h} \right)} + L^{2} + {{RL}\left( {\frac{1}{\tan \left( {j + i + h} \right)} - \frac{1}{\tan \; h}} \right)}} = 0$ ${n \times \left\{ {L - \frac{R}{\tan \; h} + \frac{R}{\tan \left( {j + i + h} \right)} + \frac{{L\; \tan \; j\; \tan \; i} + {R\left( {{\tan \; i} + {\tan \; j}} \right)}}{\tan \; j\; {\tan \left( {j + i + h} \right)}}} \right\}} + L^{2} + {{RL}\left\{ \frac{{\tan \; h} - {\tan \left( {j + i + h} \right)}}{{\tan \left( {j + i + h} \right)}\tan \; h} \right\}}$ ${{n\left\{ {L - \frac{R}{\tan \; h} + \frac{R}{\tan \left( {j + i + h} \right)} + \frac{{L\; \tan \; j\; \tan \; i} + {R\left( {{\tan \; i} + {\tan \; j}} \right)}}{\tan \; j\; {\tan \left( {j + i + h} \right)}}} \right\}} + {\frac{{\tan \; h} - {\tan \; \left( {j + i + h} \right)}}{\tan \; h \times {\tan \left( {j + i + h} \right)}}{RL}} + L^{2}} = 0$ ${n\; \frac{\begin{matrix} {{L\; \tan \; h\; \tan \; j\left\{ {{\tan \; \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {R\left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\; \tan \; h\; \tan \; j} - {\tan \; h\; \tan \; i}} \right\}} \end{matrix}}{\tan \; h\; \tan \; j\; {\tan \left( {j + i + h} \right)}}} = {- \frac{{{RL}\left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} + {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}}{\tan \; h\; \tan \; \left( {j + i + h} \right)}}$ $n = {- \; \frac{\tan \; {j\;\left\lbrack {{{RL} \times \left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} + {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}} \right\rbrack}}{\begin{matrix} {{L\; \tan \; h\mspace{11mu} \tan \; j \times \left\{ {{\tan \; \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {R \times \left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\; \tan \; h\; \tan \; i} - {\tan \; h\; \tan \; j}} \right\}} \end{matrix}}}$

Meanwhile, from equation (12),

${m = {\left( {L + \frac{R}{\tan \; j} + \frac{R}{\tan \; i}} \right) \times \frac{\tan \; i}{L}n}},$

thus,

$m = {{\left( {L + \frac{R}{\tan \; i} + \frac{R}{\tan \; i}} \right) \times \frac{\tan \; i}{L} \times n} = {\left( {L + \frac{R}{\tan \; i} + \frac{R}{\tan \; i}} \right) \times \frac{\tan \; i}{L} \times \frac{{- \tan}\; {j\left\lbrack {{{RL} \times \left\{ {{\tan \; h} - {\tan \; \left\{ {j + i + h} \right)}} \right\}} - {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}} \right\rbrack}}{\begin{matrix} {{L\; \tan \; h\; \tan \; i \times \left\{ {{\tan \; \left( {i + j + h} \right)} + {\tan \; i}} \right\}} -} \\ {R \times \left\{ {{\tan \; \left( {i + j + h} \right)\tan \; i} - {2\; \tan \mspace{11mu} h\; \tan \; i} - {\tan \; h\; \tan \; i}} \right\}} \end{matrix}}}}$ $m = {\left( {L + \frac{R}{{\tan \; j}\;} + \frac{R}{\tan \; i}} \right) \times \frac{\tan \; i}{L} \times \frac{{- \tan}\; {j\left\lbrack {{{RL} \times \left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} - {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}} \right\rbrack}}{\begin{matrix} {{L\; \tan \; h\; \tan \; j \times \left\{ {{\tan \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {R \times \left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\; \tan \; h\; \tan \; i} - {\tan \; h\; \tan \; j}} \right\}} \end{matrix}}}$ $m = {\left( \frac{{L\; \tan \; j\; \tan \; i} + {R\; \tan \; i} + {R\; \tan \; j}}{\tan \; j\; \tan \; i} \right) \times \frac{\tan \; i}{L} \times \frac{{- \tan}\; {j\left\lbrack {{{RL} \times \left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} + {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}} \right\rbrack}}{\begin{matrix} {{L\; \tan \; h\; \tan \; j \times \left\{ {{\tan \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {R \times \left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\; \tan \; h\; \tan \; i} - {\tan \; h\; \tan}} \right\}} \end{matrix}}}$ $m = {\left( \frac{{L\; \tan \; j\; \tan \; i} + {R\; \tan \; i} + {R\; \tan \; j}}{\tan \; j \times L} \right) \times \frac{{- \tan}\; {j\left\lbrack {{{RL} \times \left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} + {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}}} \right\rbrack}}{\begin{matrix} {{L\; \tan \; h\; \tan \; j \times \left\{ {{\tan \; \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {R \times \left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\tan \; h\; \tan \; i} - {\tan \; h\; \tan \; j}} \right\}} \end{matrix}}}$ $\mspace{20mu} {m = {- \frac{\left( {{L\; \tan \; j\; \tan \; i} + {R\; \tan \; i} + {R\; \tan \; j}} \right)\begin{bmatrix} {{{RL} \times \left\{ {{\tan \; h} - {\tan \left( {j + i + h} \right)}} \right\}} +} \\ {L^{2}\tan \; h\; {\tan \left( {j + i + h} \right)}} \end{bmatrix}}{\begin{matrix} {{L^{2}\tan \; h\; \tan \; j \times \left\{ \; {{\tan \left( {j + i + h} \right)} + {\tan \; i}} \right\}} -} \\ {{LR} \times \left\{ {{{\tan \left( {j + i + h} \right)}\tan \; j} - {2\; \tan \; h\; \tan \; i} - {\tan \; h\; \tan \; j}} \right\}} \end{matrix}}}}$

That is, m and n may be expressed by using known constants L (distance between rings) and R (diameter of ring), and tan h, tan i, tan j, tan (j+i+h).

Meanwhile, ∠GPH=h, ∠HPI=i, and ∠IPJ=j may be obtained from the lengths of respective sides of triangles ΔGPH, ΔHPI, and ΔHPJ.

Summarize using the second raw of cosines:

${\cos \; A} = {\frac{b^{2} + c^{2} - a^{2}}{2{bc}}.}$

Then,

${\cos \; h} = \frac{\left( \overset{\_}{PG} \right)^{2} + \left( \overset{\_}{PH} \right)^{2} - \left( \overset{\_}{GH} \right)^{2}}{2\overset{\_}{PG} \times \overset{\_}{PH}}$ ${\cos \; i} = \frac{\left( \overset{\_}{PH} \right)^{2} + \left( \overset{\_}{PI} \right)^{2} - \left( \overset{\_}{HI} \right)^{2}}{2\overset{\_}{PH} \times \overset{\_}{PI}}$ ${\cos \; j} = {\frac{\left( \overset{\_}{PI} \right)^{2} + \left( \overset{\_}{PJ} \right)^{2} - \left( \overset{\_}{IJ} \right)^{2}}{2\overset{\_}{PI} \times \overset{\_}{PJ}}.}$

Meanwhile, since the coordinates of points P, G, H, I, and J are known, the distance of each straight line may be obtained using the formula between two points.

Through these, h, i, and j values are obtained and tangent values thereof may be obtained.

Through these, m and n may be obtained.

As described above, three-dimensional spatial coordinates of points G′, H′, I′, and J′ may be obtained, and points O_(a) and O_(b) which bisect line segments G′H′ and I′J′ are calculated to obtain the second reality three-dimensional spatial coordinates of the image of the medical device.

In addition, three-dimensional spatial coordinates of the points on the ring, which are indicated by markers, may be calculated using three-dimensional spatial coordinates of points G′, H′, I′, J′, O_(a), and O_(b).

FIG. 11 shows the medical device including a marker indicating one point on the ring, the X-ray light source P, and an image of the medical device on the detector, and here, the center point of the ring is designated as Oa, the intersection point of the ring and the marker is designated as M′, and the image thereof is designated as M. Points G, H, I, and J of FIG. 10 are shown.

Point M′ is on straight line PM passing through points P and M, and distanced from the center point O_(a) of the circle by R/2 since point M′ is on the ring.

Assume that P (p₁, p₂, p₃) and M (m₁, m₂, m₃) and O_(a) (a₁, a₂, a₃) and O_(b) (b₁, b₂, b₃). Then, the equation of PM is given by

$\frac{x - p_{1}}{m_{1} - p_{1\;}} = {\frac{y - p_{2}}{m_{2} - p_{2\;}} = {\frac{z - p_{3}}{m_{3} - p_{3}}.}}$

Assume that m₁−p₁=α, m₂−p₂=β, and m₃−p₃=γ, and coordinates of M′ are (m′₁, m′₂, m′₃). Then,

$\begin{matrix} {\frac{m_{1}^{\prime} - p_{1}}{\alpha} = {\frac{m_{2}^{\prime} - p_{2}}{\beta} = {\frac{m_{3}^{\prime} - p_{3}}{\gamma}->}}} & {{Equation}\mspace{14mu} (13)} \end{matrix}$

The distance between M′ and O_(a) is R/2, and thus,

$\begin{matrix} {{\left( {b_{1} - m_{1}^{\prime}} \right)^{2} + \left( {b_{2} - m_{2}^{\prime}} \right)^{2} + \left( {b_{3} - m_{3}^{\prime}} \right)^{2}} = {\frac{R^{2}}{4}->}} & {{Equation}\mspace{14mu} (14)} \end{matrix}$

In addition, since point M′ is one point on the circle which has point O_(a) as the center point and passes through points A′, B′, G′, and H′, point M′ and the circle are present on the same plane. The line connecting points O _(a) and O_(b) is a normal line to this plane.

The equation of this plane is given by

${\frac{x - a_{1}}{a_{1} - b_{1\;}} + \frac{y - a_{2}}{a_{2} - b_{2\;}} + \frac{z - a_{3}}{a_{3} - b_{3\;}}} = 0$

Therefore,

$\begin{matrix} {{\frac{m_{1}^{\prime} - a_{1}}{a_{1} - b_{1\;}} + \frac{m_{2}^{\prime} - a_{2}}{a_{2} - b_{2\;}} + \frac{m_{3}^{\prime} - a_{3}}{a_{3} - b_{3}}} = {0->}} & {{Equation}\mspace{14mu} (15)} \end{matrix}$

By summarizing equations (13), (14), and (15),

M′ (m′₁, m′₂, m′₃) can be obtained.

∠M′O_(a)G′ and ∠M′O_(a)H′ can be obtained from the obtained coordinates of M′, G′, O_(a), and H′, and are used to find a relative position of the marker, thereby easily operating the medical device.

In cases where the rings are parallel with each other without a flexure, the position of the medical device between the rings can be found, and thus a three-dimensional image can be implemented.

In addition, in cases where a marker capable of indicating a specific site on the ring is used, the relative position of an internal portion of the medical device can be found.

At least one marker may be attached on the ring as needed. In cases where the medical device may have a flexure portion, two rings as above are attached before and after the flexure portion, thereby obtaining a three-dimensional flexible image. The number of flexure portions may be two or more.

In cases where one light source instead of a plurality of light sources is used in X-ray fluoroscopy, a predetermined marker is needed on the medical device. In addition, even in cases where a plurality of light sources are used, the use of a predetermined marker may be a supplement help in calculating the reality three-dimensional spatial coordinates of the medical device.

The predetermined marker is an identifier capable of recognizing the rotation, direction, and the like, of the medical device. In cases where a plurality of light sources are used as described above, the rotation, direction, and the like, of the medical device can be recognized. However, in cases where the three-dimensional spatial coordinates of the medical device are obtained by using one light source, it is difficult to find information about accurate rotation or direction of the medical device.

In cases where one light source is used, X-rays are irradiated to the medical device using an X-ray light source, and then an X-ray fluoroscopy image of the medical device through the X-ray light source is detected by an X-ray detector. After that, the three-dimensional spatial coordinates of the medical device are obtained based on the image of the medical device, which is detected by the X-ray detector.

The medical device of the present invention may include a radio wave reflector.

In this case, a first microwave generator generates a microwave toward the radio wave reflector.

Here, the microwave needs to have a frequency band which is harmless to humans or various medical devices in a surgical field.

In addition, a second microwave generator generates a microwave toward the radio wave reflector. The microwave by the second microwave generator and the microwave by the first microwave generator have different frequency bands, and should not cause mutual interference. Preferably, the microwave by the second microwave generator and the microwave by the first microwave generator are mutually orthogonal to each other.

After that, a microwave receiver receives a first microwave and a second microwave, which are reflected through the radio wave reflector, separates a radio wave by the first microwave and a radio wave by the second microwave through filtering, and measuring the time during which the microwaves are transmitted, reflected, and received, thereby calculating the position of the radio wave reflector, that is, the position of the medical device.

Step (c): Determining Third Reality Three-Dimensional Spatial Coordinates of Periphery of Medical Device

A three-dimensional image of the periphery of the medical device is obtained from the imaging unit of the medical device, and third reality three-dimensional spatial coordinates thereof are determined.

According to an embodiment of the present invention, the third reality three-dimensional spatial coordinates are on an orthogonal coordinate system.

The three-dimensional image of the periphery of the medical device may be obtained using a plurality of (e.g., two) cameras.

According to an embodiment of the present invention, the camera as the imaging unit includes two rings and markers attached on one-sided portions of the respective rings, as an identification part. The rings may be positioned on the same plane or a different plane with respect to an image sensor of the camera. Even though the rings are positioned on a different plane from the image sensor, the three-dimensional position of the image sensor can be obtained as long as the image sensor is positioned on the plane parallel with a plane including the rings and the distance between the two planes is known. The markers may be attached to the rings so as to coincide with the directions of the x-axis and the y-axis of the image sensor.

As described above, when the identification part is attached to the imaging unit of the medical device, the reality three-dimensional coordinates of the three-dimensional image obtained by the imaging unit can be calculated by using the identification part. In addition, the imaging unit may be positioned in vivo or ex vivo.

First, a three-dimensional endoscope using a 3D camera is exemplified. The 3D camera is composed of two cameras. When a pair of rings and markers indicating points on the rings are displayed on each of the cameras in directions of x-axis and y-axis of the image sensor, the center axis of the lens of each of the cameras and x- and y-axial vectors of the image sensor can be calculated by X-ray fluoroscopy, and based on these, the reality three-dimensional spatial coordinates of the three-dimensional image of the 3D camera can be calculated.

More specifically, the respective cameras are designated as camera A (image sensor A of FIG. 13) and camera B (image sensor B of FIG. 13). Assume that the straight line constituting the center axis of the lens of camera A is represented by Z_(a), the straight line constituting the horizontal axis on the plane of the image sensor of camera A is represented by X_(a), the straight line constituting the vertical axis is represented by Y_(a), and the intersection point of X_(a) and Y_(a) is represented by O_(a). Then, straight line Z_(a), is perpendicular to X_(a) and Z_(a) while passing through O_(a). Assume that, when the image of subject M is captured, the image (point) of the subject captured by camera A is represented by M_(a)′, the one-to-one match point on the image sensor is represented by M_(a) and the focal point of camera A is represented by F_(a). Then, M is the point on the line passing through F_(a) and M_(a).

The center axis of the lens of the camera, as described above, can be calculated from radiation images of a pair of rings installed on an external surface of the cylinder-shaped camera, and the vectors of X_(a) and Y_(a) can be calculated using the markers indicating points on the rings. Furthermore, the coordinates of the center point of the image sensor, that is, O_(a), can be calculated. When the vectors of X_(a) and Y_(a) are known and the reality three-dimensional spatial coordinates of O_(a) are known, the reality three-dimensional spatial coordinates of M_(a) can be calculated by the method illustrated in FIG. 1. In addition, the reality three-dimensional spatial coordinates of F_(a) are positioned on the center axis of the lens, and can be mechanically obtained. The equation of straight line M_(a)F_(a) can be calculated by using the reality three-dimensional spatial coordinates of points F_(a) and M_(a). Since point M is positioned on the extension line of the straight line, the equation of straight line M_(a)F_(a) is identical to the equation of straight line M_(a)M.

Since coordinates of O_(b), M_(b), and F_(b) of camera B can be found in the same manner, the equation of the straight line passing through points M and M_(b) can be obtained. In addition, since point M is the intersection point of the straight line passing through points M and M_(a) and the straight line passing through points M and M_(b), the coordinates of point M can be obtained.

Therefore, the image of the 3D camera can be made into reality three-dimensional spatial coordinates, which are then linked with the reality three-dimensional spatial coordinates of the different three-dimensional image, which has been previously obtained by X-ray fluoroscopy, thereby realizing a visual image.

As shown in FIG. 12, straight line connecting the centers of two rings is allowed to coincide with the center axis of the lens, and while specific sites of the ring, indicated by the markers, are on the same plane as the image sensor (camera A or B), the points are connected to the center point of the image sensor, thereby showing vectors of the x-axis and the y-axis of the image sensor. In addition, the markers corresponding to the x-axis and the y-axis are differently displayed for differentiation. For example, one maker is displayed longer and the other marker is displayed shorter.

Through the method described above, the center axis of the two rings is calculated and three-dimensional coordinates of the marker are calculated, thereby finding the vectors of the x axis and y axis of the image sensor and the center axis (z axis) of the lens and the center point of the image sensor.

Referring to FIG. 13, if the focal point of camera A is F_(a) at the shooting, the image M_(a) of subject M is the intersection point of the straight line passing through M and F_(a) and the image sensor.

Therefore, M is on the straight line passing through F_(a) and M_(a). For the same reason, M is on the straight line passing through F_(b) and M_(b).

In FIG. 14, three-dimensional coordinates of the center points of the image sensors, O_(a) and O_(b), can be calculated from radiation images of two pairs of rings installed on external surfaces of the respective cameras (described in FIG. 10).

In addition, the markers on the ring are displayed on the x-axis and the y-axis of the image sensor, and the coordinates of the specific sites on the ring, indicated by the markers, can be calculated by the radiation images, and thus the three-dimensional vectors of the x-axis and the y-axis of the image of each camera (described in FIG. 11).

In addition, the focal point of each camera is on the line which is perpendicular to the x-axis and the y-axis of the image sensor and passes through the center point O_(a) or O_(b) of the image sensor, that is, the line passing through the center points of one pair of rings, and the focal distance is determined by controlling the lens, and thus the three-dimensional coordinates of focal point F_(a) or F_(b) of each camera can be obtained.

Meanwhile, subject M is the intersection point of the straight line passing through M_(a) and F_(a) and the straight line passing through M_(b) and F_(b), and the three-dimensional coordinates of points M_(a), F_(b), M_(b), and F_(b) are known, and thus the point at which the two lines intersect in the three-dimensional space can be obtained.

More specifically, as shown in FIG. 14,

Assume that M (m₁, m₂, m₃), M_(a) (a₁, a₂, a₃), M_(b) (b₁, b₂, b₃), F_(a) (c₁, c₂, c₃), and F_(b) (d₁, d₂, d₃),

a₁−c₁=α₁, a₂−c₂=α₂, a₃−c₃=α₃, and

b₁−d₁=β₁, b₂−d₂=β₂, b₃−d₃−β₃

Then, straight line MM_(a) is

${\frac{x - a_{1}}{\alpha_{1\;}} = {\frac{y - a_{2}}{\alpha_{2}} = \frac{z - a_{3}}{\alpha_{3\;}}}},$

and

straight line MM_(b) is

$\frac{x - b_{1}}{\beta_{1\;}} = {\frac{y - b_{2}}{\beta_{2\;}} = {\frac{z - b_{3}}{\beta_{3}}.}}$

Since M (m₁, m₂, m₃) is the intersection point of the two straight lines,

$\frac{m_{1} - a_{1}}{\alpha_{1}} = {\frac{m_{2} - a_{2}}{\alpha_{2\;}} = \frac{m_{3} - a_{3}}{\alpha_{3\;}}}$ $m_{2} = {{\frac{\alpha_{2}}{\alpha_{1\;}} \times \left( {m_{1} - a_{1}} \right)} + a_{2}}$ $\frac{m_{1} - b_{1}}{\beta_{1}} = {\frac{m_{2} - b_{2}}{\beta_{2\;}} = \frac{m_{3} - b_{3}}{\beta_{3\;}}}$ $m_{2} = {{\frac{\beta_{2}}{\beta_{1}} \times \left( {m_{1} - b_{2}} \right)} + b_{2}}$ ${{\frac{\alpha_{2}}{\alpha_{1\;}} \times \left( {m_{1} - a_{1}} \right)} + a_{2}} = {{{\frac{\beta_{2}}{\beta_{1}} \times \left( {m_{1} - b_{1}} \right)} + {{b_{2}\left( {\frac{\alpha_{2}}{\alpha_{1\;}} - \frac{\beta_{2}}{\beta_{1}}} \right)}m_{1}}} = {{\frac{\alpha_{2}}{\alpha_{1\;}}a_{1}} - {\frac{\beta_{2}}{\beta_{1\;}}b_{1}} - a_{2} + b_{2}}}$ $m_{1} = \frac{{\frac{\alpha_{2}}{\alpha_{1\;}}a_{1}} - {\frac{\beta_{2}}{\beta_{1\;}}b_{1}} - a_{2} + b_{2}}{\frac{\alpha_{2}}{\alpha_{1}} - \frac{\beta_{2}}{\beta_{1\;}}}$ $m_{1} = \frac{\frac{{\alpha_{2}\beta_{1}a_{1}} + {\alpha_{1}\beta_{2}b_{1}} - {\alpha_{1}\beta_{1}a_{2}} + {\alpha_{1}\beta_{1}b_{2}}}{\alpha_{1}\beta_{1}}}{\frac{{\alpha_{2}\beta_{1}} + {\alpha_{1}\beta_{2\;}}}{\alpha_{1}\beta_{1}}}$ $m_{1} = \frac{{\alpha_{2}\beta_{1}a_{1}} + {\alpha_{1}\beta_{2\;}b_{2}} - {\alpha_{1}\beta_{1}a_{2}} + {\alpha_{1}\beta_{1}b_{2}}}{{\alpha_{2}\beta_{1}} + {\alpha_{1}\beta_{2\;}}}$

Here, m₂ and m₃ can be obtained by putting these values in the equation of straight line.

That is, the reality three-dimensional coordinates of the three-dimensional image, M (m₁, m₂, m₃) of the 3D camera can be obtained, and thus the three-dimensional image on the three-dimensional coordinates can be implemented.

Take an ultrasonography equipment as an example of the imaging unit. When an ultrasonic probe is radiographed, if the position of the scan line of the radiographed ultrasonic probe is made into three-dimensional spatial coordinates, the three-dimensional spatial coordinates of this image can be made into reality three-dimensional spatial coordinates since the ultrasonography image is based on the scan line.

Referring to FIG. 16, rings, which have the same center axis and are parallel with each other, are mounted on both ends of the ultrasonic probe, and markers displaying a scan line are formed on outer edges of the rings, so that the reality three-dimensional spatial coordinates of the scan line can be obtained by the method as described above. The ultrasonography image expresses a relative position from the scan line, and thus, based on this, the three-dimensional coordinates of the ultrasonography image can be made into reality three-dimensional spatial coordinates. If the ultrasonography image has three dimensions, the ultrasonography image can be three-dimensionally expressed, and if the ultrasonography image has two dimensions, the ultrasonography image can be two-dimensionally expressed in the three-dimensional image. This image may be implemented in combination with another three-dimensional image.

Step (d): Examining Images on Same Coordinates in Three Three-Dimensional Spatial Coordinates

As described above, the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body, the second reality three-dimensional spatial coordinates of the image of the medical device, and the third reality three-dimensional spatial coordinates of the three-dimensional image of the periphery of the medical device are determined, and then images on the same coordinates in the three three-dimensional spatial coordinates are examined.

Step (e): Producing Complex Reality Three-Dimensional Image

If one image is on the same coordinates in the three three-dimensional spatial coordinates, the image is selected, and if a plurality of images are on the same coordinates, necessary images among the plurality of images are selected to produce a complex reality three-dimensional image. Thus, an accurate procedure can be performed in real time using an effect of a kind of augmented reality.

According to an embodiment of the present invention, the method of the present invention further includes a step of determining a three-dimensional image of a periodic-mobility human body part (e.g., heart), which has periodic mobility and generates signals (e.g., electrocardiogram, blood pressure, organ pressure, and abdominal pressure) linked with the periodic mobility, and fourth reality spatial coordinates thereof. Step (d) is performed by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, the third reality three-dimensional spatial coordinates, and the fourth reality spatial coordinates. For example, the surface and the internal structure of the heart are uniformly changed according to the period of the heart. In addition, the change is uniform according to the position of the waveform of the electrocardiogram. Therefore, the morphological structure of the heart is uniform according to the position of the waveform of the electrocardiogram. Thus, the three-dimensional image of the heart according to the position of the waveform of the electrocardiogram is obtained (using, for example, esophagus and heart ultrasonic waves), and the three-dimensional coordinates of the image are calculated when a patient is placed on the operation table, and as a result, the real-time three-dimensional image of the heart can be implemented in linkage with the electrocardiogram even without obtaining the three-dimensional image of the heart in real time. According to another embodiment, a complex image of the heart is configured by obtaining a three-dimensional image of CT or MRI scans according to the electrocardiogram outside the operating place, obtaining a two- or three-dimensional image of heart ultrasonic waves according to the electrocardiogram, and matching these images using the same region based on the coordinates. Here, there are several complex images of the heart depending on the period of the electrocardiogram. After that, the reality three-dimensional spatial coordinates of the heart using esophageal echocardiography are allowed to coincide with the three-dimensional spatial coordinates of the complex image of the heart on the operation table, to obtain fourth real spatial coordinates of the complex image of the heart, thereby finding fourth real spatial coordinates of the complex image of the heart according to the period of the electrocardiogram. The complex reality three-dimensional image may be implemented using the esophageal ultrasonography image and the period of the electrocardiogram, or separately using the esophageal ultrasonography image and the period of the electrocardiogram.

System for Producing Complex Reality Three-Dimensional Image

A system of the present invention includes a first processor for calculating (i) first reality three-dimensional spatial coordinates of a three-dimensional image of the human body obtained before a medical device is applied to a subject; (ii) second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject; and (iii) third reality three-dimensional spatial coordinates of a three-dimensional image of the periphery of the medical device, by using the same coordinate system. Further, the system of the present invention includes a second processor for examining images on the same coordinates in the three three-dimensional spatial coordinates, by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial, and the third reality three-dimensional spatial, and then, if one image is on the same coordinates, selecting the one image, and if a plurality of images are on the same coordinates, necessary image(s) are selected from the plurality of images, thereby producing a complex reality three-dimensional image.

The system of the present invention is divided into the first processor and the second processor as elements, for definitely describing the present invention, but substantially, the two processors may be manufactured as one processor. Therefore, it would be obvious to a person skilled in the art that a processor capable of performing all of the above processes is also included in the scope of the present invention.

The system of the present invention may further include a computer-readable storage medium coupled with the processors. The storage medium is not particularly limited, and thus includes various storage media known in the art, for example, CD-R, CD-ROM, DVD, data signal embedded in a carrier wave, flash memory, floppy disk, hard drive, magnetic tape, MINIDISC, nonvolatile memory card, EEPROM, optical disk, optical storage medium, RAM, ROM, system memory, and a Web server, but is not limited thereto.

The system of the present invention may be constructed in various manners. For example, the system of the present invention may be constructed of a multi-processor computer array, a Web server, a multi-user/interactive system, or the like.

The system of the present invention may include various elements, and for example, the system of the present invention may constructed to include a data base storing information about images and spatial coordinates, a processor calculating reality three-dimensional spatial coordinates using the same coordinate system, and a processor selecting necessary image(s) to produce a complex reality three-dimensional image. 

1. A method for producing a complex reality three-dimensional image, the method comprising: (a) determining first reality three-dimensional spatial coordinates of a three-dimensional image of the human body obtained before a medical device is applied to a subject; (b) determining second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject, the medical device including an imaging unit; (c) obtaining a three-dimensional image of the periphery of the medical device from the imaging unit of the medical device, and determining third reality three-dimensional spatial coordinates of the obtained image, the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates being determined by using the same coordinate system; (d) examining images on the same coordinates in the three three-dimensional spatial coordinates by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates; and (e) producing a complex reality three-dimensional image by selecting, if one image is on the same coordinates, the one image, and selecting, if multiple images are on the same coordinates, necessary image(s) from the multiple images.
 2. The method of claim 1, wherein the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates are on an orthogonal coordinate system.
 3. The method of claim 1, wherein the same coordinate system used to determine the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates is set based on an operation table on which the medical device is applied to the subject.
 4. The method of claim 1, wherein the three-dimensional image of the human body used in step (a) is an X-ray fluoroscopy image, a computed tomography (CT) image, a magnetic resonance imaging (MRI) image, a positron emission tomography (PET) image, a PET/CT image, a PET/MRI image, a radioisotope imaging (RI) image, an ultrasonography image, or a complex image thereof.
 5. The method of claim 1, wherein the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates are determined using electromagnetic radiation, the electromagnetic radiation being X-ray radiation or radiofrequency electromagnetic radiation.
 6. The method of claim 5, wherein the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates are determined using X-ray fluoroscopy.
 7. The method of claim 5, wherein the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates are determined using mono-plan fluoroscopy, bi-plan fluoroscopy, or multi-plan fluoroscopy.
 8. The method of claim 1, wherein the first reality three-dimensional spatial coordinates of step (a) are determined using three-dimensional spatial coordinates of the three-dimensional image of the human body obtained before the medical device is applied to the subject and reality three-dimensional spatial coordinates of a particular part of the subject human body.
 9. The method of claim 8, wherein step (a) comprises the following sub-steps of: (a-1) irradiating X-rays toward the particular part of the human body of the subject using a first X-ray light source; (a-2) detecting an X-ray fluoroscopy image of the human body by the first X-ray light source, using a first X-ray detector; (a-3) irradiating X-rays toward the particular part of the human body to cross the X-rays by the first X-ray light source, using a second X-ray light source; (a-4) detecting an X-ray fluoroscopy image of the human body by the second X-ray light source, using a second X-ray detector; and (a-5) determining the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body by using the human body images detected by the first X-ray detector and the second X-ray detector.
 10. The method of claim 9, wherein sub-step (a-5) comprises: (a-5-1) determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector, by using physical three-dimensional coordinate information of the first X-ray light source and the first X-ray detector, distance information therebetween, and size information of the particular part of the human body detected by the first X-ray detector, and determining three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector, by using physical three-dimensional coordinate information of the second X-ray light source and the second X-ray detector, distance information therebetween, and size information of the particular portion of the human body detected by the second X-ray detector; (a-5-2) determining three-dimensional spatial coordinates of the particular part of the human body, by using the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the first X-ray light source and the first X-ray detector and the three-dimensional spatial coordinates of the three-dimensional image of the particular part of the human body positioned between the second X-ray light source and the second X-ray detector; and (a-5-3) allowing the three-dimensional spatial coordinates of the particular part of the human body to coincide with the three-dimensional spatial coordinates of the three-dimensional image of the human body of step (a) to determine the first reality three-dimensional spatial coordinates of the three-dimensional image of the human body.
 11. The method of claim 1, wherein the medical device includes an identification part for recognizing the position of the medical device.
 12. The method of claim 11, wherein the identification part includes at least one of: (i) at least one sphere-shaped identification part, (ii) at least one ring-shaped identification part, (iii) an identification part including a ring and at least one marker indicating a specific site of the ring, (iv) an identification part including at least one pair of rings, (v) an identification part including at least one pair of rings and at least one marker indicating a specific site of the rings, or (vi) an identification part including markers respectively indicating three differentiated points, the three points of (vi) the identification part being differentiated by using markers with different images, indicating the respective points, or indicating the points using ends of linear markers having differentiated images or intersection points of the markers.
 13. The method of claim 12, wherein the identification part includes two or more rings and one or more makers, the two or more rings being mounted on the medical device at different positions, the one or more markers being attached to one side of each of the rings.
 14. The method of claim 13, wherein the two or more rings have different imaging characteristics.
 15. The method of claim 13, wherein in cases where the number of markers is two or more, the markers have different imaging characteristics.
 16. The method of claim 12, wherein the second reality three-dimensional spatial coordinates are determined using one or two or more X-ray fluoroscopy images, and wherein step (b) comprises the following sub-steps of: (b-1) irradiating X-rays toward the medical device using an X-ray light source; (b-2) detecting an X-ray fluoroscopy image of the medical device by the X-ray light source, by using an X-ray detector; and (b-3) determining the second reality three-dimensional spatial coordinates of the image of the medical device, by using the through-position of the identification part, which is seen through on the X-ray fluoroscopy image.
 17. The method of claim 16, wherein in sub-step (b-3), the second reality three-dimensional spatial coordinates of the image of the medical device and three-dimensional spatial coordinates of the marker are determined by using the through-position of the identification part, which is seen through on the X-ray fluoroscopy image, thereby grasping the relative position of an internal portion of the medical device from the three-dimensional spatial coordinates of the marker.
 18. The method of claim 16, wherein sub-step (b-3) is performed using three-dimensional position information of the identification part, or three-dimensional position information and dimension information of the identification part.
 19. The method of claim 1, wherein the medical device includes a radio wave reflector.
 20. The method of claim 19, wherein step (b) comprises the following sub-steps of: (b-1) generating a first microwave toward the radio wave reflector using a first microwave generator; (b-2) generating a second microwave toward the radio wave reflector using a second microwave generator; (b-3), after sub-steps (b-1) and (b-2), receiving the first and second microwaves reflected from the radio wave reflector; and (b-4) determining the second reality three-dimensional spatial coordinates of the image of the medical device by using information of the received first and second microwaves.
 21. The method of claim 1, wherein the imaging unit includes a plurality of cameras, the three-dimensional image of the periphery of the medical device being obtained using the plurality of cameras in step (c).
 22. The method of claim 21, wherein the camera as the imaging unit includes an identification part.
 23. The method of claim 1, wherein the imaging unit is an ultrasonic probe, the three-dimensional image of the periphery of the medical device being obtained using the ultrasonic probe in step (c).
 24. The method of claim 23, wherein the ultrasonic probe as the imaging unit includes an identification part.
 25. The method of claim 1, wherein the three-dimensional image of the human body of step (a) and the three-dimensional image of the periphery of the medical device of step (c) are different kinds of images.
 26. The method of claim 1, further comprising a step of determining a three-dimensional image of a periodic-mobility human body part, which has periodic mobility and generates signals linked with the periodic mobility, and fourth reality spatial coordinates thereof, wherein step (d) is performed by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, the third reality three-dimensional spatial coordinates, and the fourth reality spatial coordinates.
 27. A system for producing a complex reality three-dimensional image, the system comprising: (a) a first process for calculating first reality three-dimensional spatial coordinates of a three-dimensional image of the human body obtained before a medical device is applied to a subject; (ii) second reality three-dimensional spatial coordinates of an image of the medical device applied to the subject; and (iii) third reality three-dimensional spatial coordinates of a three-dimensional image of the periphery of the medical device, by using the same coordinate system; and (b) a second processor for examining images on the same coordinates in the three three-dimensional spatial coordinates, by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial, and the third reality three-dimensional spatial, and then, selecting, if one image is on the same coordinates, the one image, and selecting, if a plurality of images are on the same coordinates, necessary image(s) from the plurality of images, to produce a complex reality three-dimensional image.
 28. The system of claim 27, wherein the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates are on an orthogonal coordinate system.
 29. The system of claim 27, wherein the same coordinate system used to determine the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, and the third reality three-dimensional spatial coordinates is set based on an operation table on which the medical device is applied to the subject.
 30. The system of claim 27, further comprising an X-ray fluoroscopy imaging device, a computed tomography (CT) imaging device, a magnetic resonance imaging (MRI) imaging device, a positron emission tomography (PET) imaging device, a PET/CT imaging device, a PET/MRI imaging device, a radioisotope imaging (RI) imaging device, or an ultrasonography imaging device, for obtaining the three-dimensional image of the human body.
 31. The system of claim 27, further comprising an X-ray fluoroscopy imaging device or a radiofrequency electromagnetic radiation device for obtaining the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates.
 32. The system of claim 27, wherein the first reality three-dimensional spatial coordinates and the second reality three-dimensional spatial coordinates are calculated by using information obtained from mono-plan fluoroscopy, bi-plan fluoroscopy, or multi-plan fluoroscopy.
 33. The system of claim 27, wherein the first reality three-dimensional spatial coordinates are determined using three-dimensional spatial coordinates of the three-dimensional image of the human body obtained before the medical device is applied to the subject, and reality three-dimensional spatial coordinates of a particular part of the subject human body.
 34. The system of claim 33, wherein the first reality three-dimensional spatial coordinates are calculated using information obtained from mono- or bi-plan fluoroscopy, the information including: (i) physical three-dimensional coordinate information of a first X-ray light source and a first X-ray detector, distance information therebetween, and size information of the particular part of the human body detected by the first X-ray detector; (ii) physical three-dimensional coordinate information of a second X-ray light source and a second X-ray detector, distance information therebetween, and size information of the particular part of the human body detected by the second X-ray detector; and (iii) three-dimensional spatial coordinates of the three-dimensional image of the human body.
 35. The system of claim 27, further comprising a medical device, wherein the medical device includes an identification part for recognizing the position of the medical device.
 36. The system of claim 35, wherein the identification part includes at least one of: (i) at least one sphere-shaped identification part, (ii) at least one ring-shaped identification part, (iii) an identification part including a ring and at least one marker indicating a specific site of the ring, (iv) an identification part including at least one pair of rings, (v) an identification part including at least one pair of rings and at least one marker indicating specific sites of the rings, or (vi) an identification part including markers respectively indicating three differentiated points, the three points of (vi) the identification part being differentiated by using markers with different images, indicating the respective points, or indicating the points using ends of linear markers having differentiated images or intersection points of the markers.
 37. The system of claim 36, wherein the identification part includes two or more rings and one or more makers, the two or more rings being mounted on the medical device at different positions, the one or more markers being attached to one side of each of the rings.
 38. The system of claim 37, wherein the two or more rings have different imaging characteristics.
 39. The system of claim 37, wherein in cases where the number of markers is two or more, the markers have different imaging characteristics.
 40. The system of claim 36, wherein the second reality three-dimensional spatial coordinates are calculated by using the through-positions of the rings of the identification part, which are seen through on an X-ray fluoroscopy image.
 41. The system of claim 40, wherein the second reality three-dimensional spatial coordinates are calculated by using the through-positions of the rings and markers of the identification part, which are seen through on the X-ray fluoroscopy image, and thus the second reality three-dimensional spatial coordinates of the image of the medical device and three-dimensional spatial coordinates of the marker are determined, thereby obtaining the relative position of an internal portion of the medical device from the three-dimensional spatial coordinates of the markers.
 42. The system of claim 41, wherein the second reality three-dimensional spatial coordinates are calculated by using position information of the identification part, or position information and dimension information of the identification part.
 43. The system of claim 35, wherein the medical device includes a radio wave reflector.
 44. The system of claim 35, wherein the medical device includes an imaging unit.
 45. The system of claim 44, wherein the imaging unit includes a plurality of cameras, the three-dimensional image of the periphery of the medical device being obtained using the plurality of cameras.
 46. The system of claim 44, wherein the camera as the imaging unit includes an identification part.
 47. The system of claim 44, wherein the imaging unit is an ultrasonic probe.
 48. The system of claim 47, wherein the ultrasonic probe as the imaging unit includes an identification part.
 49. The system of claim 27, wherein the three-dimensional image of the human body and the three-dimensional image of the periphery of the medical device are different kinds of images.
 50. The system of claim 27, further comprising a processor for calculating a three-dimensional image of a periodic-mobility human body part, which has periodic mobility and generates signals linked with the periodic mobility, and fourth reality spatial coordinates thereof, wherein the second processor produces a complex reality three-dimensional image by using the first reality three-dimensional spatial coordinates, the second reality three-dimensional spatial coordinates, the third reality three-dimensional spatial coordinates, and the fourth reality spatial coordinates. 